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Abstract 


This  report  describes  the  research  performed  by  the  Colorado  Center  for 
Astrodynamics  Research  (CCAR)  under  an  AASERT  Grant  to  study  Global 
Positioning  System  based  attitude  determination  for  spinning  satellites. 
Two  key  contributions  are  discussed  -  phase-based  methods  for  estimation 
of  the  spin  axis  of  a  spinning  satellite,  and  signal-to-noise  ratio  (SNR) 
based  methods  for  estimation  of  one  or  more  spacecraft  axes.  The  latter 
method,  which  is  very  robust  and  widely  applicable,  is  demonstrated  on 
actual  flight  data  from  three  NASA  missions.  It  provides  instantaneous 
pointing  knowledge  at  the  3-10  deg  level  with  a  single  GPS  antenna  and  is 
useful  for  attitude  initialization,  verification,  GPS  phase  ambiguity  resolu¬ 
tion,  and  a  reliable  backup  attitude  determination  sensor. 
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CHAPTER  1 


INTRODUCTION 

This  report  describes  the  research  performed  by  the  Colorado  Cen¬ 
ter  for  Astrodynamics  Research  (CCAR)  under  an  AASERT  Grant  to  study 
Global  Positioning  System  based  attitude  determination  for  spinning  satel¬ 
lites.  This  research  augments  a  series  of  parent  grants  from  the  Naval 
Research  Laboratory  on  GPS-based  attitude  estimation.  The  subsequent 
chapters  of  this  report  are  extracted  from  the  doctoral  dissertation  of 
Charles  Behre  submitted  to  the  Graduate  School  of  the  University  of  Colo¬ 
rado,  Boulder,  in  December  1997,  entitled  “GPS  Based  Attitude  Algorithms 
for  Low  Cost  Satellite  Missions.” 

1.1  Research  Motivation 

The  original  goal  of  this  research  was  to  develop  algorithms  for 
using  GPS  onboard  spin  stabilized  satellites  ranging  from  low-cost  student 
projects  such  as  the  Student  Nitric  Oxide  Explorer  (SNOE)  to  highly 
maneuverable,  advanced  Navy  satellites  such  as  Living  Plume  Shield 
(LIPS)  III.  In  the  process  of  investigating  traditional  phase-based  tech¬ 
niques  we  came  across  suggestions  for  the  use  of  an  alternate  GPS  observ¬ 
able,  the  signal-to-noise  ratio  for  attitude  estimation.  Our  research  has 
fully  developed  this  type  of  approach  which  has  very  broad  applicability  to 
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not  only  spacecraft  attitude  determination  but  for  land  vehicle  orientation 
as  well. 

1.2  Satellite  Attitude  Determination 

The  GPS  system  consists  of  a  constellation  of  24  Earth  orbiting  sat¬ 
ellites  that  transmit  ranging  signals  and  a  ground  control  facility.  Recep¬ 
tion  of  the  signals  provides  a  variety  of  information.  Typically,  receivers 
track  quantities  such  as  the  ranges  to  the  satellites  in  view,  the  carrier 
phases  of  the  signals,  and  the  signal-to-noise  ratios  (SNR)  of  the  measure¬ 
ments.  While  GPS  is  most  recognized  as  a  tool  for  position  estimation,  it  is 
also  widely  used  for  attitude  estimation. 

Many  types  of  instruments  are  used  for  satellite  attitude  determina¬ 
tion.  They  include  star  sensors,  sun  sensors,  Earth  horizon  sensors,  gyros, 
magnetometers,  and  GPS  receivers.  To  minimize  cost,  the  choice  of  atti¬ 
tude  sensor  has  to  be  well  matched  with  the  mission  requirements.  Choos¬ 
ing  a  sensor  with  a  much  higher  level  of  performance  than  is  required 
wastes  resources.  Star  sensors  have  accuracies  of  better  than  0.1  deg,  but 
they  are  heavier,  more  expensive,  and  require  more  power  than  other  atti¬ 
tude  sensors  [Wertz,  1978].  For  accuracies  of  under  1  deg,  companies  such 
as  Ithaco  and  Barnes  produce  fine  sun  sensors  and  Earth  horizon  sensors. 
These  instruments  are  priced  in  the  range  of  fifty  to  over  a  hundred  thou¬ 
sand  dollars.  Coarse  magnetometers  provide  accuracies  of  approximately 
3  to  10  degrees.  The  prices  range  from  ten  to  twenty  thousand  dollars 
depending  upon  the  desired  accuracy  and  the  number  of  axis  to  be  sensed. 
Coarse  sun  sensors  are  usually  the  cheapest  of  the  instruments.  When 
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bought  commercially,  they  cost  about  five  thousand  dollars,  however,  they 
must  be  used  in  conjunction  with  another  sensor  [Wertz  and  Larson,  1991]. 

Several  student-built  satellites,  that  have  yet  to  be  launch,  utilize 
various  combinations  of  lower  cost  instruments  to  provide  attitude  estima¬ 
tion  of  about  5  degrees.  The  spin  stabilized  SNOE  satellite,  built  at  the 
University  of  Colorado,  uses  two  horizon  crossing  indicators  and  a  three- 
axis  magnetometer  [Barth  et  al.,  1994].  The  three-axis  stabilized  CATSAT 
from  the  University  of  New  Hampshire  has  a  fiber  optic  gyro  (FOG),  a 
three-axis  magnetometer,  and  a  GPS  receiver  [CATSAT  Home  Page].  The 
gravity  gradient  stabilized  ASUSAT  satellite  from  the  Arizona  State  Uni¬ 
versity  uses  an  array  of  student  built  Earth  and  Sun  sensors  [ASUSAT  1]. 
This  last  satellite  is  of  particular  interest  to  this  research  because  it  also 
has  a  GPS  receiver  on-board  for  navigation  information. 

1.3  Previous  Work 

The  most  common  form  of  GPS-based  attitude  estimation  involves 
the  utilization  of  specially  designed  receivers  that  use  carrier  phase  inter¬ 
ferometry  between  several  GPS  antennas.  Performance  of  1  deg  or  better 
has  been  shown  using  GPS  attitude  receivers  such  as  the  Trimble  Vector 
[Cohen,  1992],  [Trimble,  1993]  and  the  Ashtech  3DF  [Lu,  et  al.,  1993],  as 
well  a  system  of  separate  receiver  modules  [Lu,  et  al.,  1993],  [Cannon,  et 
al.,  1994]. 

The  use  of  GPS  for  spacecraft  attitude  determination  onboard  an 
Earth-pointing  spacecraft  already  has  been  demonstrated  on  the  U.S.  Air 
Force  RADCAL  (Radar  Calibration)  satellite  [Lightsey  et  al.,  1994],  [Ward 
and  Axelrad,  1995],  on  the  Crista  Spas  mission  [Ward,  L.M.,  1996],  [Brock 
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et  al.,  1995],  on  the  REXII  satellite  [Lightsey  et  al.,  1996],  and  on  the  space 
shuttle  GANE  (GPS  Attitude  and  Navigation  Experiment)  mission  [Car¬ 
penter  and  Hain,  1997].  Each  of  these  missions  use  a  Trimble  TANS  Vec¬ 
tor  receiver.  The  Vector  is  capable  of  tracking  up  to  six  GPS  satellites  on 
four  antennas  and  measuring  carrier  phase  differences. 

Martin-Niera  and  Lucas  of  the  European  Space  Agency  were  the 
first  to  suggest  the  adaptation  of  GPS  phase  measurements  specifically  to 
spinning  satellites  [Martin-Niera  and  Lucas,  1992].  They  describe  a  fast 
Fourier  transform  (FFT)  based  method  for  using  a  single  GPS  antenna 
baseline  to  determine  spacecraft  spin  and  nutation  rates  as  well  as  nuta¬ 
tion  angles  and  orientation  of  the  angular  momentum  vector.  Their 
approach  uses  triple-differenced  phase  measurements  as  the  basic  observ¬ 
able.  Their  experimental  results  are  based  on  data  collected  by  a  ground 
GPS  receiver  mounted  on  a  spinning  platform.  The  spin  frequencies  cho¬ 
sen  for  the  simulation  minimize  the  sampling  errors  inherent  in  the  FFT 
technique. 

The  receivers  that  are  able  to  measure  differential  phase  tend  to  be 
more  expensive  than  single  antenna  navigation  receivers.  However,  atti¬ 
tude  estimation  from  a  single  antenna  is  also  possible  by  exploiting 
another  measured  quantity.  The  SNR  measurement,  that  is  commonly 
computed  by  all  GPS  receivers  as  they  track  incoming  signals,  is  also 
dependent  upon  vehicle  attitude.  Due  to  the  variations  in  the  GPS 
antenna  gain  pattern,  the  SNR  values  are  related  to  the  orientation  of  the 
antenna  with  respect  to  the  GPS  satellites  which  are  transmitting  the  par¬ 
ticular  signals  being  received.  In  this  research,  the  nature  of  these  varia¬ 
tions  is  exploited  to  estimate  the  boresight  of  the  antenna. 
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Single  axis  estimation  using  GPS  SNR  data  was  first  suggested  in  a 
Motorola  patent  [Ghaem  et  al.,  1992].  The  inventors  propose  combining  a 
directional  antenna  with  a  hemispherical  antenna  in  a  GPS  receiver  to 
determine  heading.  The  GPS  receiver  indicates  the  direction  to  at  least 
one  GPS  satellite.  The  directional  antenna  is  then  rotated  until  it  receives 
a  signal  from  the  chosen  satellite.  The  combination  of  the  orientation  of 
the  directional  antenna  and  the  known  coordinates  of  the  GPS  satellite 
establish  a  reference  heading.  The  user’s  heading  is  computed  from  the 
angle  between  the  directional  antenna  and  the  travel  direction.  No  perfor¬ 
mance  estimates  are  given  in  the  patent. 

In  another  patent,  Flawn  suggests  using  a  multi-element  antenna  in 
conjunction  with  a  GPS  receiver  [Flawn,  1996] .  From  a  common  GPS  sat¬ 
ellite,  the  various  antenna  elements  have  a  different  reception  of  the 
incoming  signal.  This  difference  is  then  utilized  to  determine  the  orienta¬ 
tion  of  the  antenna  with  respect  to  the  GPS  satellite  and  as  a  result,  the 
users  heading. 

A  proposal  for  using  SNR  measurements  from  a  standard  antenna 
for  spacecraft  attitude  determination  is  made  by  Hashida  and  Unwin  of 
the  University  of  Surrey.  They  suggest  using  the  SNR  measurement  from 
one  antenna  for  directional  information  [Hashida  and  Unwin,  1993]. 
These  readings  compose  one  set  of  the  observations  in  a  Kalman  filter.  The 
filter  estimates  position  and  velocity  from  GPS  data  as  well  as  attitude 
from  the  SNR  data  of  a  single  antenna.  Their  approach  is  to  model  both 
the  azimuth  and  elevation  pattern  of  an  antenna  from  ground  data  and 
then  to  apply  this  model  in  the  filter.  Accuracies  of  3  deg  in  roll  and  pitch 
and  10  deg  in  yaw  are  shown  only  on  simulated  data. 
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Another  European  group  uses  the  SNR  readings  from  two  canted 
antennas  to  estimate  the  body  frame  coordinates  of  the  line-of-sight  (LOS) 
vectors  from  the  GPS  antenna  to  the  GPS  satellites  [Serrano,  et  al.,  1995]. 
The  relationship  between  the  SNR  and  the  boresight  vector  of  the  antenna 
is  modeled  in  detail.  Their  results  are  demonstrated  entirely  with  simu¬ 
lated  data  and  depended  greatly  on  the  degree  of  SNR  mismodeling. 

Researchers  at  the  Jet  Propulsion  Laboratory  use  the  SNR  measure¬ 
ments  as  weights  in  a  LOS  averaging  method  for  coarse  attitude  estima¬ 
tion  [Dunn,  1997] .  They  estimate  the  boresight  vector  of  a  GPS  antenna  by 
averaging  all  of  the  LOS  vectors  to  the  visible  GPS  satellites.  The  LOS 
vectors  are  weighted  such  that  the  ones  corresponding  to  the  largest  SNR’s 
are  assigned  the  highest  weights.  The  estimate  is  used  on  the  GPS/MET 
satellite  for  the  purpose  of  satellite  acquisition. 

1.4  Research  Contribution 

This  research  extends  the  work  on  GPS  based  attitude  determina¬ 
tion  in  two  key  areas.  The  first  is  the  use  of  GPS  phase  data  for  spinning 
satellites.  The  second  is  a  very  robust  and  widely  applicable  technique 
based  on  SNR  measurements  made  from  a  single  antenna.  The  former 
algorithms  are  demonstrated  for  simulated  data  and  a  ground  experiment 
conducted  with  the  assistance  of  the  Naval  Research  Laboratory.  The  lat¬ 
ter  algorithms  are  demonstrated  using  actual  flight  data  from  three  NASA 
satellite  missions  as  well  as  ground  test  data. 

The  phase-based  approach  includes  several  algorithms.  For  finding 
the  angular  rates,  a  simple  averaging  technique,  FFT  technique,  and  an 
auto-regressive  (AR)  method  are  studied.  For  determining  the  orientation 
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of  the  angular  momentum  vector,  an  averaging  method  and  a  batch  filter 
are  analyzed.  Results  are  shown  on  a  computer  simulation  of  GPS  phase 
measurements  collected  on  a  spinning  satellite  and  on  data  collected  from 
a  ground  test.  The  computer  simulation  models  a  perfectly  symmetrical 
spinning  satellite  with  two  GPS  antennas  mounted  along  a  diameter  of  the 
face.  The  ground  experiment  consists  of  a  tipping  and  spinning  platform 
that  approximates  the  kinematics  of  a  spinning  satellite.  Data  is  taken 
from  a  Trimble  TANS  Vector  receiver. 

The  SNR-based  approach  uses  a  maximum  likelihood  estimation 
(MLE)  of  an  antenna  boresight  vector  given  a  set  of  estimated  off-boresight 
angles.  Results  are  shown  on  a  variety  of  data  sets.  These  sets  include 
space  based  data  and  ground  based  data  with  different  GPS  receivers  and 
antennas.  Extension  of  the  single  axis  based  estimation  to  both  spinning 
and  three-axis  stabilized  satellites  is  also  examined. 

1.5  Research  Overview 

The  second  chapter  of  this  report  presents  a  background  of  GPS 
phase  based  measurements  and  their  relationship  to  attitude.  This  is  fol¬ 
lowed  by  a  detailed  description  of  the  use  of  carrier  phase  measurements 
applied  specifically  to  spinning  satellites. 

The  remaining  chapters  deal  specifically  with  the  SNR-based  esti¬ 
mation.  Chapter  4  describes  the  relationship  between  the  SNR  measure¬ 
ment  and  the  antenna  boresight  vector.  The  first  part  of  Chapter  5 
presents  an  estimation  method  for  a  single  spacecraft  axis  using  the  maxi¬ 
mum  likelihood  estimate  (MLE)  of  an  antenna  boresight  vector.  The  sec¬ 
ond  part  examines  the  use  of  this  vector  for  various  scenarios  of  attitude 
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determination.  Chapter  6  details  the  various  data  sets  that  are  used  to 
analyze  the  performance  of  the  MLE  technique.  Chapter  7  presents  and 
examines  the  results  of  the  MLE  approach  applied  to  these  sets. 


CHAPTER  2 


GPS  BASED  ATTITUDE  BACKGROUND 

2.1  Introduction 

This  chapter  describes  the  phase  measurement  made  by  GPS  receiv¬ 
ers  and  how  these  measurements  are  used  to  compute  vehicle  attitude. 
The  first  part  details  the  computation  of  single  difference  phase  measure¬ 
ments  made  between  two  antennas  and  time  differenced  phase  measure¬ 
ments.  These  quantities  are  related  to  the  baseline  vector  between  the 
antennas  and  to  the  displacement  of  the  baseline  vector  between  two 
epochs,  respectively.  The  next  section  presents  the  equations  and  quanti¬ 
ties  necessary  to  estimate  the  baselines  and  the  baseline  displacements. 
Finally,  the  baselines  are  related  to  three-axis  and  single-axis  vehicle  atti¬ 
tude. 

2.2  Phase  Difference  Measurements 

The  GPS  satellites  transmit  an  encoded  carrier  signal  at 
1575.42  MHz  corresponding  to  a  wavelength  of  19.03  cm.  To  track  this 
incoming  signal,  a  GPS  receiver,  using  the  code  correlation  technique, 
must  generate  a  similarly  encoded  reference  signal  at  the  same  frequency 
[Hoffman- Wellenhof,  et.  al.,  1992],  In  this  method,  the  reference  signal  is 
shifted  with  respect  to  time  so  that  the  two  signals  optimally  match.  The 
fractional  phase  difference  is  computed  by 
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<|>  =  -f5t,  (2.1) 

where  <J>  is  the  difference  between  the  carrier  phases  of  the  two  signals,  f  is 
the  carrier  frequency,  and  5t  is  the  time  between  signal  transmission  and 
reception.  This  is  known  as  the  carrier  beat  phase.  The  unknown  number 
of  cycle  shifts  needed  for  perfect  cycle  synchronization  is  called  the  integer 
ambiguity.  As  the  GPS  satellites  and  receiving  antennas  move,  the 
required  time  or  Doppler  shift  changes.  To  maintain  lock  on  the  signal  the 
receiver  must  continuously  adjust  the  phase  shift  of  the  reference  signal. 
As  long  as  the  receiver  remains  locked,  the  change  in  phase  is  tracked  over 
any  number  of  additional  complete  cycles  or  reduction  in  complete  cycles. 
This  is  known  as  the  “integrated  Doppler”  or  “accumulated  delta  phase” 
[Axelrad  and  Brown,  1996],  However,  if  the  receiver  momentarily  loses 
lock  on  the  signal,  the  full  cycle  count  becomes  invalid  and  therefore,  track¬ 
ing  must  be  restarted  at  the  current  phase  fraction.  This  loss  of  signal  lock 
is  referred  to  as  a  cycle  slip. 

By  using  two  antennas,  the  difference  between  the  phase  of  the 
arriving  carrier  signal  between  the  two  antennas  can  be  measured,  as 
denoted  by  a  A  in  this  research.  As  before,  the  measured  difference  is  only 
a  fractional  measurement.  The  number  of  complete  cycles  between  the  two 
antennas  is  unknown.  If  the  antennas  are  operating  under  a  common 
receiver  clock,  then  the  phase  difference  between  the  two  antennas  is 
related  to  the  baseline  vector  connecting  the  antennas  by 

A<\>  =  (rL)TeL-k  +  p  +  v  (2.2) 

where  r  is  the  baseline  vector,  e  is  the  line  of  sight  to  the  GPS  satellite,  k 
is  the  integer  ambiguity,  p  is  the  antenna  line  bias  in  fractions  of  a  cycle, 
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and  v  is  additional  random  noise.  The  line  bias  represents  the  difference 
in  the  delays  between  the  two  antennas  caused  by  cable  and  electronic  dif¬ 
ferences.  The  superscript  L  refers  to  the  local  system  in  which  the  line-of- 
sight  (LOS)  vectors  are  computed.  Figure  2.1  show  this  relationship 
graphically. 


GPS  antennas 


Figure  2.1.  Relationship  between  an  antenna  phase  difference  and  a 
baseline  vector. 


A  third  type  of  phase  measurement  used  in  this  research  is  the  time 
differenced  measurement,  denoted  by  a  8.  It  is  formed  by  differencing  two 
A<|>  measurements  taken  with  the  same  antenna  pair  at  two  distinct 
epochs.  The  time  differenced  phase  measurement  between  epoch  t,  and  tk 
is  defined  as 
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8A<|>(ti)sA<Kti  +  k)-A<|>(ti).  (2.3) 

Another  type  of  phase  difference  that  is  not  used  in  this  research  is 
between  measurements  from  two  different  satellites.  This  is  commonly 
denoted  with  the  symbol  V . 

2.3  Baseline  Estimation  from  A<p  Measurements 

Given  two  or  more  baseline  vectors  on  a  vehicle  provided  by  three  or 
more  antennas,  the  attitude  of  the  vehicle  can  be  computed  using  GPS 
phase  measurements.  The  baseline  vectors  are  computed  from  a  batch  of 
A<J>  measurements  taken  from  multiple  GPS  satellites.  After  neglecting  the 
noise,  M  phase  measurements  from  a  GPS  satellite  i  are  related  to  the  M 
baselines  by 

=  e^B  +  kj,  (2.4) 

where 

KiA4>i2...A*iM],  (2.5) 

B  *  [rl  r2  •••  rM]  •  (2-6) 

and 

=  [(-kil  +  Pl)  (—  ^i2  **”  P2)  •••  +  (2.7) 

The  integer  ambiguities  are  found  from  the  integer  parts  of  the  elements  of 
k  and  the  line  biases  are  found  from  the  common  fractional  parts.  Each 
GPS  satellite  and  baseline  pair  has  a  unique  integer  ambiguity  associated 
with  it.  However,  the  line  bias  are  unique  only  to  each  baseline.  If  mea- 
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surements  from  N  GPS  satellites  are  taken  simultaneously,  then  an  equa¬ 
tion  relating  the  phase  differences  to  the  baselines,  LOS  vectors,  integers, 
and  line  biases  is 


A<D  =  H 


where 


(2.8) 


T 

A<D 


(2.9) 


e2  ...  dN 
1  0  ...  0 
0  1  : 

:  \  0 
0  ...  0  1 


(2.10) 


and 


,T,T 
kl  k2 


(2.11) 


The  baselines,  integers,  and  biases,  are  found  by  solving  the  normal  equa¬ 
tion  in  (2.8)  for  the  matrices  B  and  K. 

In  the  matrix  B  every  baseline  corresponds  to  three  unknown  vector 
components.  Every  GPS  satellite  creates  one  unknown  integer  ambiguity 
per  baseline  in  the  matrix  K  and  provides  one  phase  measurement  per 
baseline.  This  leaves  more  unknowns  than  equations.  However,  if  the 
GPS  receivers  maintain  a  continuous  track  on  the  satellites,  the  integer 
ambiguities  remain  fixed  over  a  number  of  measurements  with  only  the 
baseline  components  changing.  If  there  is  a  sufficient  enough  change  in 
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the  H  matrix,  then  this  quickly  leads  to  more  measurements  than 
unknowns,  making  Equation  (2.8)  overdetermined.  For  E  measurement 
epochs,  M  baselines,  and  N  GPS  satellites,  the  requirement  for  Equation 
(2.8)  to  be  overdetermined  is 

(3ME  +  MN)  <  MNE .  (2.12) 

For  the  ratio  of  unknowns  to  measurements  to  be  less  than  or  equal  to  one 
the  condition  is 


3E+N 

NE 


<1. 


(2.13) 


When  N=l,2,3  or  E=1  the  ratio  in  Equation  (2.13)  is  never  smaller  than 
one  for  any  value  of  E.  Therefore,  a  direct  solution  for  both  the  baselines 
and  integers  requires  a  minimum  of  four  satellites  and  two  measurement 
epochs. 

Much  larger  data  sets  are  needed  for  reliable  estimates.  For  exam¬ 
ple,  to  accurately  determine  the  baselines  with  the  TANS  Vector,  the  user’s 
manual  recommends  collecting  data  for  a  duration  of  at  least  8  hours  taken 
every  2  to  5  seconds  [Trimble,  1993].  This  is  primarily  to  allow  enough 
time  for  a  sufficient  change  in  the  LOS  vector  components  of  the  H  matrix. 
Furthermore,  cycle  slips  and  the  inclusion  of  measurements  from  GPS  sat¬ 
ellites  acquired  in  the  middle  of  the  data  set  create  additional  integer 
ambiguities  that  have  to  be  computed. 


2.4  Displacement  Vector  Estimation  from  8Aq>  Measurements 

Time  differenced  measurements  can  be  used  to  estimate  the  dis¬ 
placement  of  a  baseline  between  two  measurement  epochs,  as  illustrated 
in  Figure  2.2.  Given  two  A(j)  measurements  at  times  ti  and  t^,  the  time  dif- 
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Figure  2.2.  The  displacement  vector. 

ferenced  phase,  SA^),  can  be  computed  using  the  definition  presented  in 
Equation  (2.3).  Substituting  Equation  (2.2)  into  Equation  (2.3),  assuming 
no  cycle  slips,  and  approximating  the  line  of  sight  vector  as  constant  over 
this  interval,  an  equation  involving  the  antenna  displacement  vector,  dr , 
can  be  formed  as 

5A<J) (tj)  =e  (t/drtfj)  +  (ki  +  k-kj)  +  (vi  +  k-vi)  .  (2.14) 

When  the  receiver  remains  locked  on  to  the  satellite  signal,  the  two  integer 
ambiguities,  ki+k  and  ki;  are  equal,  resulting  in 

8A<j>(ti)  =  6  (tj)Tdr(tj)  +  (vi  +  k-vi)  .  (2.15) 

To  solve  Equation  (2.15)  for  the  three  elements  of  dr ,  measurements 
from  at  least  three  GPS  satellites  must  be  available.  Assuming  that  the 
noise  is  uncorrelated,  unbiased,  and  equi-variant,  then  the  least  squares 
solution  is 
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dr(t.)  =  (gTg)_1GT 


8A<j>i  (tj) 
8A<()2  (tj) 


|5A<()n  (t.) 


(2.16) 


where 

G=ks2...6N]T,  (217) 

the  subscript  refers  to  each  satellite,  and  N  is  the  number  of  satellites.  If  a 
cycle  slip  occurs  between  the  measurement  of  A<])(ti)  and  A^t^),  then  the 
integer  ambiguities  at  the  two  epochs  are  not  equal  and  Equation  (2.14)  is 
not  solvable.  However,  each  cycle  slip  only  affects  the  computation  of  a  sin¬ 
gle  displacement  vector,  making  it  a  very  robust  observable. 

In  general,  the  use  of  this  type  of  measurement  usually  does  not  pro¬ 
duce  high  accuracy  estimates  because  the  measurement  noise  is  higher 
and  the  geometry  is  weaker  than  the  A(f>  or  VA<{>  observable.  As  a  result, 
they  are  normally  used  for  initialization  such  as  in  surveying  to  establish 
an  initial  estimate  for  long  antenna  baselines  [Hoffman- Wellenhof,  et.  al., 
1992]  and  in  attitude  determination  as  a  means  for  initial  ambiguity  reso¬ 
lution  [Cohen,  1996] .  However,  if  there  is  a  large  and  rapid  movement  of 
the  baseline,  then  these  two  conditions  are  no  longer  a  factor. 

2.5  Attitude  Estimation  from  Vector  Observations 

B  L 

The  attitude  of  the  vehicle  is  defined  as  the  matrix,  C  ,  that 
relates  a  vector  in  a  given  local  frame  denoted  by  the  superscript  L  to  a 
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representation  in  a  vehicle  body  fixed  frame  denoted  by  the  superscript  B. 
The  representations  of  the  vector,  r,  in  the  two  frames  are  related  by 

rB  =  BCLrL,  (2.18) 

where  BCL  is  the  product  of  three  Euler  rotations  computed  by 


BCL  =  Rx  (6)  R2  (v)  R3  (*)  • 


(2.19) 


Given  two  non-parallel  vectors,  ri  and  r2,  known  in  the  vehicle  fixed 
frame  and  measured  in  the  local  frame,  an  attitude  matrix  is  determin¬ 
able.  Using  the  Gram-Schmidt  process  [Fraleigh  and  Beauregard,  1990]  a 
set  of  orthonormal  vectors,  ej,  e2,  and  e3  is  computed  by  the  equations 


ei  = 


1  » 


e2  " 


(  L 

L 

r2  *  el 

r9  - 

ei 

Z 

^el  *  el ) 

(2.20) 


and 


e3  -  ei  x  e2 


The  resulting  local  to  body  attitude  matrix  is 


BCL  = 


(2.21) 


For  more  than  two  observations,  the  noise  in  the  measurements 
allow  for  multiple  attitude  matrices  to  be  computed  depending  upon  the 
choice  of  and  r2.  The  best  estimate  of  the  attitude,  minimizes  the 
Wahba  cost  function 
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L(BCL)  =  \  X  ak|rB-BCLrL|,  (2.22) 

k=  1 

where  n  is  the  number  of  non-parallel  vector  observations  and  ak  are  the 
measurement  weights.  The  solution  to  Equation  (2.22)  is  the  orthonormal 
matrix  that  best  fits  the  measurements  in  a  weighted  least  squares  sense. 
Many  solutions  to  Equation  (2.22)  have  been  previously  offered  [Markley, 
1993],  [Shuster,  1981],  [Wahba,  1965]. 

2.6  Attitude  Estimation  from  Phase  Measurements 

After  collecting  a  batch  of  differential  phase  measurements  over  sev¬ 
eral  epochs,  the  baselines  in  the  local  frame  and  the  integer  ambiguities 
are  computed  by  solving  Equation  (2.8).  Assuming  the  GPS  receiver  stays 
continuously  locked  on  to  the  GPS  satellite,  the  integer  ambiguities  only 
have  to  be  determined  once.  New  ambiguities  are  readily  computed  from 
the  current  information.  After  solving  for  the  integers,  the  baselines  are 
more  easily  estimated  as  they  change.  Rearrange  Equation  (2.8)  to  form 

A<D-K  =  HB  (2.23) 

where  the  H  matrix  is  simplified  to 

HT=  [e,  e2...  e„],  (2.24) 


Solving  for  the  baselines  produces 
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B  =  (HTH)  1HT(AO-K)  .  (2.25) 

In  Equation  (2.25)  there  are  3M  unknowns  at  every  measurement  epoch. 
Every  GPS  satellite  gives  M  observations.  Therefore,  at  least  three  satel¬ 
lites  must  be  tracked  to  solve  for  the  baselines  explicitly  at  each  epoch. 

If  at  least  two  baselines  are  known  in  the  body  fixed  frame,  the 
results  of  Equation  (2.25)  are  used  to  estimate  the  vehicle  attitude.  For 
two  baselines,  this  attitude  is  computed  by  Equations  (2.20)  and  (2.21). 
Multiple  baselines  require  the  minimization  of  Equation  (2.22).  Figure  2.3 
shows  an  example  representation  of  the  antenna  baselines  in  both  the  local 
and  body  frames.  In  this  case  the  vehicle  body  y-axis  is  defined  along  one 
of  the  antenna  baselines  and  the  xy-body  plane,  the  shaded  region  in  the 
figure,  is  set  as  the  plane  containing  both  baselines. 


Figure  2.3.  Relationship  between  antenna  baselines  and  reference 
frames. 
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2.7  Single-Axis  Attitude 

As  shown  in  Equation  (2.19),  three  angles  are  required  to  completely 
define  the  vehicle  attitude  and  thus,  the  three  orthonormal  axes.  Single¬ 
axis  attitude  requires  just  two  distinct  angles.  These  two  angles  define  the 
orientation  of  one  of  the  body  axes  with  respect  to  the  reference  frame  as 
shown  in  Figure  2.4.  The  plane  containing  the  other  two  axes  is  known  to 
be  perpendicular  to  the  first  axis.  However,  the  orientation  of  the  remain¬ 
ing  two  axes  within  the  plane  is  not  specified.  Single-axis  attitude  is  useful 
when  the  pointing  direction  of  a  single  instrument  is  needed. 


Figure  2.4.  Single-axis  attitude. 

One  example  of  a  vehicle  with  the  requirement  for  only  single-axis 
attitude  knowledge  is  a  spinning  satellite.  As  described  in  more  detail  in 
Chapter  3,  a  spinning  satellite  has  only  one  inertially  fixed  axis.  The  other 
two  are  continuously  rotating  about  this  axis.  The  placement  of  GPS 
antennas  on  the  plane  containing  these  rotating  axes  creates  a  baseline 
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vector  that  undergoes  both  large  and  rapid  displacement.  As  a  result,  the 
problems  with  using  the  displacement  vector  observable,  as  discussed  at 
the  end  of  section  2.4,  are  not  a  significant  factor.  This  measurement  is  the 
basis  for  the  attitude  estimation  techniques  presented  in  the  next  chapter. 

Another  area  where  single-axis  estimation  is  useful  is  for  integer 
ambiguity  resolution.  As  discussed  in  Sections  2.3  and  2.6,  baseline  and 
attitude  estimation  using  phase  measurements  have  an  ambiguity  that 
must  determined.  This  is  accomplished  when  either  a  certain  number  of 
measurements  are  available  or  by  using  a  search  techniques  to  evaluate  all 
of  the  possible  integer  combinations.  The  performance  of  the  latter  method 
is  improved  if  the  search  space  is  narrowed.  The  phase-based  method 
described  by  Ward  [Ward,  L.M.,  1996],  requires  a  single  axis  to  be  known  to 
within  20  degrees  for  a  rapid  determination  of  the  attitude  and  integer 
ambiguities. 

2.8  Conclusions 

A  variety  of  phase  difference  measurements  can  be  computed  from 
the  GPS  signal.  The  basic  phase  measurement,  <j>,  is  formed  by  synchroniz¬ 
ing  the  GPS  receiver  reference  signal  with  the  incoming  satellite  signal. 
The  observable,  A<>,  is  made  by  differencing  the  phase  measurements 
between  two  antennas.  Differencing  this  value  between  two  epochs  forms 
8A<{). 

The  complete  vehicle  attitude  with  respect  to  a  local  system  is  speci¬ 
fied  by  determining  the  components  of  two  or  more  body  axes  in  this  sys¬ 
tem.  A  GPS  receiver  may  estimate  the  body  axes  in  the  local  system 
through  the  use  of  A<j>  measurements.  This  process  consists  of  first  deter- 
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mining  the  integer  ambiguities  of  the  carrier  phase  between  antennas, 
estimating  the  receiver  line  biases,  and  the  baseline  components.  After¬ 
wards,  only  the  baseline  components  need  to  be  computed. 

Single-axis  estimation  has  fewer  unknowns  to  estimate  at  the  cost  of 
a  less  specific  representation  of  attitude.  The  8A<|)  measurements  are  use¬ 
ful  in  determining  attitude  for  this  case.  The  next  chapter  explores  a  tech¬ 
nique  specifically  for  spinning  satellites.  The  estimation  process  is  shown 
to  be  simpler  and  have  fewer  requirements  than  for  three-axis  attitude. 
One  significant  advantage  is  that  there  is  no  need  for  ambiguity  resolution. 


CHAPTER  3 


GPS  PHASE-BASED  SPIN  AXIS  ESTIMATION 

3.1  Introduction 

In  this  chapter  several  methods  are  described  that  use  the  displace¬ 
ment  vector  observable,  as  presented  in  the  previous  chapter,  to  determine 
the  angular  rates  and  the  spin  axis  of  a  rotating  satellite.  For  angular  rate 
estimation,  two  frequency  domain  and  one  time  domain  technique  are  pre¬ 
sented.  The  first  frequency  domain  method  is  based  on  the  Fast  Fourier 
Transform  (FFT),  the  second  is  based  on  the  Auto-Regressive  (AR)  system 
model.  The  time  domain  method  involves  an  averaging  of  the  displace¬ 
ment  vectors.  For  spin  axis  estimation,  an  averaging  technique  using  dis¬ 
placement  vector  cross  products  and  a  batch  filter  are  presented.  It  is 
shown  that  the  AR  model  is  the  best  approach  for  estimating  the  satellite 
frequencies.  The  cross  product  method  is  the  simplest  technique  for  esti¬ 
mating  the  spin  axis  and  it  works  as  well  as  the  batch  filter  algorithm. 
They  both  have  estimation  errors  of  under  0.1  deg  after  400  measure¬ 
ments.  The  batch  filter  has  the  advantage  of  being  able  to  estimate  more 
parameters  than  the  cross  product  method. 

This  chapter  begins  with  a  description  of  an  ideal  spinning  satellite 
model  and  the  simulated  data  generated  from  it.  This  is  followed  by  a 
description  of  a  ground-based,  hardware  simulation  of  the  satellite  motion. 
An  actual  GPS  receiver  is  used  to  collect  the  experimental  data.  For  this 
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work,  both  the  computer  and  hardware  simulations  are  based  on  a  receiver 
that  generates  phase  difference  observations  (A<J>)  between  two  antennas 
and  computes  appropriate  line-of- sight  vectors  (e )  to  the  GPS  satellites 
based  on  its  current  navigation  solution.  Note  that  the  receiver  is  assumed 
to  be  operating  in  a  steady  state  condition  -  that  is,  it  has  already  acquired 
at  least  four  satellites  and  is  computing  position  solutions  in  real-time. 
The  acquisition  process  typically  takes  up  to  30  seconds  for  a  receiver  that 
has  been  initialized  with  an  a  priori  estimate  of  the  orbit,  GPS  ephemeris, 
and  time,  and  up  to  30  minutes  for  a  receiver  in  "cold-start".  The  remain¬ 
ing  sections  of  this  chapter  describe  the  frequency  and  attitude  estimation 
algorithms  and  the  corresponding  results  from  using  both  the  simulated 
and  experimental  data. 

3.2  Spinning  Satellite  Model 

Figure  3.1  illustrates  the  important  quantities  describing  the  kine¬ 
matics  of  a  spinning  satellite.  For  simplicity,  the  vehicle  is  assumed  to  be 
axisymetric  about  the  z®-axis.  The  master  antenna  is  located  at  the  center 
with  the  master-slave  baseline  perpendicular  to  the  axis  of  symmetry.  The 
satellite  rotates  about  the  zB-axis  at  a  rate  of  0)p.  Simultaneously,  the  zB- 
flvis  rotates  about  the  angular  momentum  axis,  H ,  at  the  rate  to,  at  an 
angle  of  0.  The  satellite  has  a  moment  of  inertia  about  the  zB-axis  of  Is  and 
moment  of  inertia  about  the  xB  and  yB-axes  of  It. 

The  body  xB-axis  is  defined  by  the  position  of  the  slave  antenna.  The 
relationships  between  to,  top,  to h  0,  and  the  inertias  are 
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instantaneous 
rotation  axis 


Figure  3.1.  Kinematic  model. 

a  =  !ilio>,cose  (3.D 

ls 

and 

to2  =  G)2  +  ft)2  +  2o)pco/cose .  (3.2) 

The  body  frame  is  fixed  to  the  rotating  satellite  and  is  denoted  with 
a  superscript  B.  The  angular  momentum  frame,  denoted  by  the  super¬ 
script  H,  is  defined  with  the  z-axis  along  the  spacecraft  angular  momen¬ 
tum  vector.  The  inertial  frame  is  designated  by  the  superscript  I.  The 
angular  momentum  vector  unit  vector,  H,  in  the  inertial  frame  is 
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i  r  iT 

[  -  [k  n  hj  ■ 


(3.3) 


3.2.1  Antenna  Motion 

In  the  absence  of  external  torques,  a  GPS  antenna  on  the  rim  of  a 
spinning  satellite  moves  according  to  the  model  described  in  Equations 
(3.1)  and  (3.2).  In  the  body  fixed  system,  the  antenna  position  remains  con¬ 
stant.  The  position  of  the  antenna  with  respect  to  the  center  of  the  space¬ 
craft  in  the  inertial  system  can  be  computed  from  the  product  of  two 
separate  rotation  matrices  given  by 


I  I~H  H„B  B 

r  =  C  C  r  , 


(3.4) 


where 


=  r0  0  0] 


(3.5) 


is  the  position  of  the  antenna  in  the  body  frame.  The  first  rotation,  C  , 
takes  r  from  the  body  to  the  angular  momentum  frame.  The  angles  y  and 
<p  are  defined  as  functions  of  the  satellite  rotation  rates  as  follows: 


V  =  o>pt  +  \|ro 

cp  =  c^t  +  cpo 


Then  HCB  can  be  expressed  as 


„  _  c\pc<p  -  c0s\j/s<p  -svirc(p  -  c0c\j/s(p  s0scp 

_ 

V'  -  cystp  +  C0S\J/C(p  -s\|/scp  +  c0cyccp  s0ccp 


s0s\p 


s0c\p 


(3.7) 


The  matrix  C  rotates  a  vector  from  the  angular  momentum  sys¬ 
tem  to  the  inertial  system.  It  is  defined  by  the  two  angles  r\  and  a  which 
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are  determined  from  the  components  of  the  angular  momentum  unit  vector 
by 


( 

q  =  cos1 


and 


a  =  cos1  (hz) . 


(3.8) 


(3.9) 


This  matrix  is 


ICH 


erica  -sa  crisa 
sr|ca  cr\  sr|sa 
— sa  0  ca 


(3.10) 


If  Equation  (3.4)  is  differenced  in  time,  an  expression  relating  the 
antenna  displacement  vectors  to  the  satellite  kinematics  can  be  formu¬ 
lated.  For  a  constant  orientation  of  a  satellite  angular  momentum  vector, 
this  relationship  at  time  t*  is 


dAt^'cVcVi.)  -HcB(ti)  i 


(3.11) 


The  displacement  vectors  in  the  local  frame  are  computed  from  the 
GPS  phase  differences  by  Equation  (2.16).  Equation  (3.11)  describes  how 
these  components  are  related  to  the  spacecraft  attitude. 

The  frequencies  observed  in  dr1  are  due  to  the  products  of  the  sines 
and  cosines  in  Equation  (3.7).  The  three  frequencies  are  (o /  ±  cop  and  cop. 
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3.2.2  Computer  Simulation 

A  computer  program  was  constructed  to  simulate  the  motion  of  a 
spinning  satellite  in  low  earth  orbit.  LOS  measurements  are  generated 
from  a  simulated  constellation  of  24  GPS  satellites  using  a  typical  set  of 
orbital  elements  and  taking  into  account  satellite  visibility  limitations. 
Perfect  phase  difference  values  are  formed  from  the  line  of  sight  and 
known  antenna  baseline  vectors.  Simulated  A<|>  measurements  are  con¬ 
structed  by  adding  zero  mean  Gaussian  noise  with  a  standard  deviation  of 
5  mm.  Displacement  vectors  are  computed  from  the  noisy  A(J)  values.  Mea¬ 
surements  are  recorded  at  2  Hz  for  12  minutes.  The  spin  rate  is  5  rpm  and 
the  nutation  angle  is  5  deg. 

In  the  remainder  of  this  chapter  the  rotation  matrices  referring  to 
the  simulated  LEO  motion  are  denoted  with  a  subscript  “s”.  Thus,  the 
expressions  given  in  Equations  (3.7)  and  (3.10)  are  for  C  and  C  , 

S  S' 

respectively. 

3.3  Ground  Experiment 

In  order  to  validate  the  results  of  the  computer  simulation,  a  series 
of  ground  tests  was  conducted.  Data  were  taken  from  a  GPS  receiver 
mounted  on  a  spin  table  located  at  a  Naval  Research  Laboratory  test  facil¬ 
ity.  The  matrices  describing  the  experimental  motion  are  denoted  with  the 
subscript  “e”. 

3.3.1  Experimental  Apparatus 

The  complete  test  apparatus  consists  of  two  pieces  of  equipment 
with  a  combined  total  of  three  rotational  axes  as  diagramed  in  Figure  3.2. 
A  photograph  of  the  actual  setup  is  shown  in  Figure  3.3.  The  TRT-7  two- 
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avis  tilter  is  mounted  on  a  support  base.  In  the  local  level  coordinate 
frame,  the  first  axis  is  oriented  along  an  east-west  direction.  The  second 
avis  is  aligned  perpendicular  to  the  first  and  rotates  around  it.  When  the 
first  rotation  angle  is  zero,  the  second  axis  lies  in  the  north-south  direction. 
The  maximum  tilt  angle  about  both  axes  is  approximately  ±45  degrees. 
The  second  piece  of  equipment  is  the  BDS-5  spinner.  It  provides  a  constant 
angular  rate  about  a  third  axis.  This  axis  rotates  according  to  the  tilt 
angles  of  the  TRT-7.  It  is  mounted  so  that  when  the  first  two  rotation 
angles  are  zero  its  orientation  is  in  the  up-down  direction.  The  GPS  anten¬ 
nas  are  mounted  on  the  top  plane  of  the  satellite  structure  perpendicular 
to  this  axis. 


Figure  3.2.  Equipment  diagram. 
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Figure  3.3.  Experimental  setup. 


3.3.2  Experimental  Motion 
H  B 

The  matrix  Cg ,  derived  in  Equation  (3.7)  is  a  3-1-3  sequence  of 
rotations  computed  from  the  Euler  angles  \jr,  cp,  and  0.  In  order  to  properly 
simulate  the  motion  of  a  spinning  satellite,  the  structure  in  the  experimen¬ 
tal  setup  has  to  be  moved  according  to  the  rotations  defined  by  the  TRT-7 
and  the  BDS-5.  For  this  case,  the  body  frame  is  defined  by  the  symmetry 
axis  of  the  structure  and  the  position  of  one  of  the  GPS  antennas  as  shown 
in  Figure  3.2.  In  the  body  frame  the  position  of  the  GPS  antenna  is  always 
on  the  xB-axis.  The  local  system  is  defined  by  the  orientation  of  the  TRT- 
7’s  tilt  axis  when  the  plane  of  the  structure  is  level  to  the  ground.  While 
there  is  no  true  angular  momentum  axis,  an  equivalent  frame  to  the  angu¬ 
lar  momentum  frame  is  defined  by  the  plane  of  the  satellite  structure  when 
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there  is  only  pure  spin.  For  the  purpose  of  comparison  to  the  simulated 
data,  this  is  referred  to  as  the  H  frame. 

If  the  east  direction  is  called  the  1-axis,  the  north  direction  the  2- 
axis,  and  the  up  direction  the  3-axis,  then  the  series  of  rotations  used  to 
move  the  GPS  antenna  is  a  1-2-3  sequence.  Defining  the  three  rotation 
angles  as  {5,  a,  and  y,  respectively,  the  total  rotation  matrix  is 

Bc"(M.y)  =B,(r)Rj(«)B1(P).  (3.12) 

In  order  to  simulate  the  correct  movement  of  a  GPS  antenna  on  a 
spinning  satellite  governed  by  Equation  (3.7),  the  angles  P,  a,  and  y,  must 
be  varied  in  time  as  a  function  of  xp,  cp,  and  0.  These  relations  are  computed 
by  setting 

Bc“  (P,  a,  Y)  =  Bcf  (fc  8,  V) .  (3.13) 

The  resulting  expressions  for  a,  P,  and  y,  are  a  function  of  the  sums 
and  products  of  the  sinusoids  of  <p,  0,  and  \|/.  These  general  equations  are 
very  involved;  however,  for  small  nutation  angles  they  may  be  approxi¬ 
mated  by 

p(t)  »0sin(<Kt)),  (3.14) 

a(t)  ~0cos  (<J>(t))  ,  (3.15) 

and 

y(t)  =<t>(t)  +V(t)  .  (3.16) 

Due  to  limitations  in  the  control  and  performance  of  the  test  equip¬ 
ment,  two  compromises  are  required.  The  first  is  to  eliminate  the  motion 
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of  p  due  to  Equation  (3.14).  The  second  is  to  approximate  the  motion  of  a 
in  Equation  (3.15)  with  a  triangle  wave.  The  equation  for  a  is 

sin  (5X(t) )  -  ...J  (3.17) 

where 


a(t)  =0-H"sin(?i(t))-isin(3X(t))  +± 

TC  L  9  25 


X(t)  =  coat  +  X0  (3.18) 

and 

co  =  —co,.  (3.19) 

a  n  1 

Combining  the  motion  of  a  about  the  2-axis  and  the  motion  of  y  about  the 
3-axis  yields  the  rotation  matrix  from  the  body  system  to  the  angular 
momentum  system: 

cycoc  -syca  sa 

Ce  =  sy  cy  0  (3.20) 

-cysa  sysa  ca 

The  resulting  motion  described  by  Equation  (3.20)  and  illustrated  in 


Figure  3.4  is  a  combination  of  spinning  and  tipping. 
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The  tipping  is  contained  in  a  single  plane  centered  on  the  He  axis. 
The  orientation  of  He  is  fixed  and  described  by  the  two  angles  o  and  r|. 
The  angle  <j  is  about  the  1-axis  and  the  angle  r|  is  about  the  2-axis.  The 
corresponding  rotation  matrix  from  the  He  frame  to  the  inertial  system  is 

\  cn  o  st, 

'-'e  =  ST|SC  C  O  -CT|Sa 
-STJCG  SO  CT|CG 

The  rotation  matrices  in  Equations  (3.20)  and  (3.21)  are  analogous 
to  the  matrices  in  Equations  (3.7)  and  (3.10)  for  the  true  satellite  rotations. 
One  difference  is  that  in  the  case  of  an  actual  spinning  satellite,  the  spin 
axis  rotates  about  the  angular  momentum  axis;  whereas  for  the  experi¬ 
mental  motion  the  spin  axis  tips  about  some  central  axis.  Another  differ¬ 
ence  is  in  the  frequency  modes  of  the  ground  experiments. 


(3.21) 
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3.3.3  Experimental  Motion  Frequencies 

For  this  case  the  two  principle  rates  are  co  and  coa  and  the  observed 
frequencies  are  at  0)  and  co±kcoa,  where  k  is  any  integer.  The  frequencies 
contained  in  the  experimental  motion  are  similar  to  those  of  a  spinning 
satellite.  It  is  shown  that  the  equations  for  satellite  motion  contain  three 
frequencies  derived  from  the  two  satellite  rates  of  cop  and  ©/.  For  the  exper¬ 
iment,  there  is  a  frequency  due  to  the  rate  about  the  spin  axis  and  frequen¬ 
cies  due  to  the  combined  rates  of  the  spin  axis  and  the  tilting  motion  about 
a  fixed  vector.  Because  a(t)  is  a  triangle  wave  it  also  contains  odd  harmon¬ 
ics  of  coa  as  opposed  to  just  a  single  frequency.  Thus  cosine  a  contains  even 
multiples  while  the  sine  a  contains  odd  ones.  A  Power  Spectral  Density 
(PSD)  plot  of  the  displacement  vectors  yields  an  infinite  number  of  spikes 
corresponding  to  all  possible  values  of  k.  However,  only  the  first  few  multi¬ 
ples  of  coa  have  enough  power  to  show  above  the  noise  floor. 

The  reference  value  for  the  ground  test  spin  rate  was  0.0850  Hz  ± 
7X10-4  Hz.  The  platform  was  leveled  to  within  1  deg  and  the  initial  align¬ 
ment  error  was  determined  to  be  ±  0.01  deg. 

3.3.4  GPS  Receiver 

The  experiments  were  conducted  using  a  Trimble  Navigation  TANS 
Vector  Attitude  Receiver  [Cohen,  1992],  [Trimble,  1993].  This  receiver  has 
6  parallel  tracking  channels  which  multiplex  among  four  antennas  to  pro¬ 
vide  attitude  estimates  as  well  as  a  navigation  solutions.  One  antenna  is 
designated  as  the  master  and  the  other  three  are  slaves.  The  receiver 
measures  the  relative  carrier  phase  for  each  slave  with  respect  to  the  mas¬ 
ter.  For  these  experiments  the  raw  differential  phase  data  is  recorded  at  a 
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rate  of  2  Hz.  Typical  differential  phase  errors  for  the  receiver  are  in  the 
range  of  2-5  mm  1-cr  depending  on  the  signal  environment  [Axelrad  and 
Chesley,  1993]. 

The  next  sections  describe  various  algorithms  for  attitude  and  fre¬ 
quency  estimation. 


3.4  Frequency  Estimation  Techniques 

In  this  section,  three  frequency  estimation  techniques  are  examined. 
The  first  method  is  a  time  domain  technique,  while  the  next  two  are  fre¬ 
quency  domain  methods. 


3.4.1  Displacement  Vector  Magnitude  Algorithm 

For  small  nutation  angles,  0)  can  be  approximated  by  averaging  the 
values  of  a  quantity  denoted  by  co*  that  is  computed  from  the  displacement 
vectors.  If  accurate  values  of  the  satellite’s  inertia  ratios  are  known,  esti¬ 
mates  of  (Dp  and  0),  can  also  be  found.  The  simplest  approximation  for  the 
angular  velocity  0)  given  in  Equation  (3.2)  can  be  computed  from  the  mag¬ 
nitude  of  dr .  Every  two  successive  positions  of  a  rotating  antenna  creates 
an  angular  displacement  centered  at  the  spin  axis.  This  angle  shown  in 
Figure  3.5  is  computed  by 


0  =  2  sin*1 


(3.22) 


where  r0  =  |r| .  Defining  co*  as  the  mean  angular  rate  of  0,  the  magnitude 
of  co*  varies  with  cop,  its  minimum  value  is 


•  I  ' 

At 


vV1* 


cop  and  its  maximum 


value  is  oi  For  a  nutation  angle  of  zero,  co=co*. 
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Figure  3.5.  Angular  displacement. 

3.4.2  FFT  Algorithm 

To  determine  the  modal  frequencies  and  modal  amplitudes  of  a  sig¬ 
nal,  the  following  steps  are  applied:  perform  a  discrete  FFT  on  the  sam¬ 
pled  signal  data;  compute  the  power  spectral  density  (PSD)  of  the 
transformed  data;  identify  the  peaks  of  the  PSD;  determine  the  frequency 
at  which  each  peak  occurs;  and  relate  these  frequencies  back  to  a  kine¬ 
matic  model. 

Using  only  the  PSD  values,  the  frequencies  are  the  most  likely  esti¬ 
mates  of  each  signal  mode.  Ideally,  they  are  the  exact  frequencies  of  the 
modes.  However,  the  sampling  properties  of  the  signal  can  lead  to  errors. 

The  frequency  resolution  of  the  DFFT  is  equal  to  the  sample  fre¬ 
quency  divided  by  the  number  of  data  points.  Ideally,  the  signal  should  be 
sampled  so  that  the  modal  frequencies  are  integer  multiples  of  the  fre¬ 
quency  bin.  In  reality,  however,  the  peaks  of  the  PSD  occur  in  the  fre¬ 
quency  bins  that  are  closest  to  the  correct  modes.  The  signal  power  of  each 
mode  is  spread  over  a  range  of  frequency  bins  surrounding  the  correct  fre¬ 
quency.  In  order  to  determine  a  better  estimate,  the  frequency  center  of 
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the  bins  is  computed  by  weighting  them  according  to  their  power.  How¬ 
ever,  the  benefit  of  this  approach  may  be  limited  by  noise. 

In  the  implementation  of  an  FFT  algorithm  there  are  no  assump¬ 
tions  about  the  structure  of  the  signal.  In  other  words,  the  FFT  can  be 
applied  to  any  type  of  signal.  There  are  other  frequency  estimation  meth¬ 
ods  which  can  take  advantage  of  the  knowledge  of  the  signal's  structure. 
The  next  section  discusses  one  of  these  parametric  methods  called  the  auto 
regressive  or  AR  method. 

3.4.3  AR  Algorithm 

In  general,  a  discrete-time  process  can  be  well  approximated  by  a 
time  series  or  rational  transfer  function  model  [Kay,  19881  •  In  particular  a 
system  of  order  2p,  can  be  represented  by  the  recursive  difference  relation 


2p 

x[n]  =  -  £  a  [k]  x  [n-k]  +  u  [n]  .  (3.23) 

k  =  1 

Taking  the  Z-transform  of  Equation  (3.23)  yields  the  transfer  function 


X(z)  _±_ 
U(z)  P(z) 


(3.24) 


where 


z 


o  + j0 
e 


(3.25) 


and 
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2p  k 

P(z)  =  1+  £  a[k]z  .  (3.26) 

k  =  i 

Given  a  sequence  of  sampled  data  x[n],  the  coefficients,  a[k],  can  be 
estimated  using  the  modified  covariance  method  [Kay,  1988].  The  roots  of 
P(z)  corresponding  to  sinusoids  are  described  by 


* 


zk  = 


±je; 

e 


(3.27) 


They  occur  in  complex  conjugate  pairs  and  lie  on  the  unit  circle  (i.e.  a  =  0) 
at  angles  corresponding  to  the  sinusoidal  frequencies  normalized  by  the 
sample  frequency: 


*  -  -  fk 

0k  =  27ifk,  fk  =  f  (3.28) 

rs 

Note  that  this  technique  directly  estimates  the  modal  frequencies  and 
therefore  does  not  provide  estimates  which  are  quantized  by  discrete  fre¬ 
quency  bins  as  does  the  FFT  method. 

The  poles  corresponding  to  the  signal  modes  lie  closest  to  the  unit 
circle,  while  the  poles  corresponding  to  the  noise  lie  farther  away.  The 
angles  of  the  signal  mode  poles  are  approximately  equal  to  27ifk .  This  is 
illustrated  in  Figure  3.6.  For  the  ideal  case  where  there  is  no  noise  the  sig¬ 
nal  is  perfectly  modeled  by  Equation  (3.26),  with  exactly  p  sinusoids,  and 
each  pole  of  Equation  (3.24)  lies  exactly  on  the  unit  circle  at  an  angle  of 
27tfk.  For  the  more  realistic  case  a  higher  order  model  than  2p  is  used  to 
account  for  the  noise  in  the  signal.  As  the  model  order  is  increased,  there 
are  more  frequencies  available  for  the  noise.  As  a  result,  the  poles  for  the 
signal  modes  get  closer  to  the  unit  circle  and  the  angles  get  closer  to  2rcfk . 
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The  model  order  is  limited  by  the  number  of  samples,  N.  A  maximum 
order  of  N/2  is  recommended  [Kay,  1988].  The  model  order  used  in  this 
work  is  70  for  values  of  N  that  are  200  and  greater. 


*:  poles  corresponding  to  signal 
x:  poles  corresponding  to  noise 

Model  order  =  10 
Number  of  sinusoids  =  2 


Figure  3.6.  The  pole  -  frequency  relationship  for  a  noisy  sampled  signal. 


3.4.4  Frequency  Estimation  Results 

The  performance  of  the  frequency  estimation  techniques  are  evalu¬ 
ated  on  both  the  simulated  and  experimental  data.  The  AR  and  displace¬ 
ment  vector  magnitude  methods  produce  frequency  errors  under  0.01  rad/s 
after  30  seconds  of  data.  The  FFT  method  can  not  resolve  the  frequencies 
with  this  small  of  a  data  set  and  still  shows  large  errors  after  100  seconds. 
For  large  batch  sizes,  the  errors  are  below  0.001  rad/s  for  both  the  AR  and 
FFT  techniques.  Averaging  ©*,  however,  does  not  show  improvement  with 
larger  amounts  of  data.  Tables  3.1  and  3.2  compare  the  performance  of  the 
three  frequency  estimation  techniques  as  a  function  of  the  data  batch  size. 
All  of  the  data  is  sampled  at  2  Hz.  For  the  simulation  data,  the  nutation 
angle  is  5  deg.  For  the  experimental  data  the  maximum  tip  angle  is  5.3 

deg.  In  the  experiment  the  programmed  value  for  coa  is  0.628  rad/s.  How- 
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ever,  due  to  the  limitations  of  the  equipment  as  noted  earlier,  this  number 
is  believed  to  be  significantly  in  error.  Judging  by  the  consistency  of  the 
methods  over  various  batch  sizes  and  by  the  simulation  results,  it  is 
believed  that  the  results  from  the  frequency  estimation  algorithms  more 
accurately  represent  the  actual  motion  than  the  programmed  values. 


Table  3.1.  Performance  of  frequency  methods  on  simulated  data. 


Method 

N  = 

100 

N  = 

200 

N  = 

800 

CO, 

(rad/s) 

C0p 

(rad/s) 

CO, 

(rad/s) 

cop 

(rad/s) 

CO, 

(rad/s) 

cop 

(rad/s) 

True  Value 

0.6789 

-0.1561 

1^1 

-0.1561 

AR 

0.6760 

-0.1530 

-0.1565 

FFT 

no  sol. 

no  sol. 

-0.1562 

Average  co* 

0.6780 

-0.1565 

BH 

0.6759 

-0.1560 

Table  3.2.  Performance  of  frequency  methods  on  experimental  data. 


N  = 

100 

N  = 

200 

N  = 

800 

Method 

CO 

(rad/s) 

coa 

(rad/s) 

CO 

(rad/s) 

(rad/s) 

CO 

(rad/s) 

(rad/s) 

Reference 

Value 

0.5349 

(1) 

0.5349 

(1) 

0.5349 

(1) 

AR 

0.5379 

0.5807 

0.5356 

0.5959 

0.5345 

0.5973 

FFT 

no  sol. 

no  sol. 

0.5381 

0.5947 

0.5350 

0.5971 

Average  co* 

0.5337 

(2) 

0.5350 

(2) 

0.5358 

(2) 

Note:  (1)  There  is  no  adequate  reference  value  for  coa. 


(2)  The  average  to*  method  does  not  produce  estimates  for  coa. 

3.5  Attitude  Estimation  Techniques 

This  section  presents  two  methods  for  determining  the  orientation  of 
the  spacecraft  angular  momentum  vector.  The  first  technique  uses  the 
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average  of  the  displacement  vector  cross  products.  The  second  uses  a 
batch  processor. 


3.5.1  Cross  Product  Algorithm 

In  the  absence  of  nutation,  a  GPS  antenna  located  on  one  face  of  a 
rotating  satellite  travels  in  a  plane  perpendicular  to  the  angular  momen¬ 
tum  axis.  The  displacement  vector  of  the  baseline  formed  by  two  moving 
antenna  is  also  in  this  plane.  Taking  the  cross  product  of  two  displacement 
vectors  yields  a  vector  with  approximately  the  same  orientation  as  the 
angular  momentum  axis.  Furthermore,  the  components  of  H  define  the 
orientation  of  the  satellite  spin  axis  in  the  inertial  frame.  Because  the 
GPS  measurements  are  taken  in  the  inertial  frame,  the  displacement  vec¬ 
tor  cross  products  are  computed  in  the  inertial  frame,  and  therefore,  H  is 

also  computed  in  the  inertial  frame. 

When  nutation  is  introduced,  the  cross  product  of  two  successive  dis¬ 
placement  vectors  is  not  the  true  orientation  of  H .  However,  as  the  satel¬ 
lite’s  instantaneous  spin  axis  rotates  about  H,  the  average  of  all  the 
vectors  formed  from  these  cross  products  tends  towards  the  actual  orienta¬ 
tion  of  H .  This  is  computed  by 


N  drCy  xdr(ti  +  k) 
H  -  N.f1|dr(ti)xdr(ti  +  k)| 


(3.29) 


where  N  is  the  number  of  measurements. 


3.5.2  Spinning  Satellite  Batch  Processor 

A  batch  processing  algorithm  for  estimation  of  the  attitude  is  also 
examined.  The  state  comprises  the  nutation  angle  0,  the  angular  position 
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of  the  antenna  about  the  body  zB-axis  vp,  the  angular  position  about  Hs ,  (p, 
and  the  two  constant  angles,  r|  and  a,  which  relate  the  inertial  to  the  angu¬ 
lar  momentum  frame.  The  resulting  state  vector  is 

X  =  [0  \j/  (p  r|  q]  •  (3.30) 

For  the  torque  free  case  considered  here,  the  state  estimate  is  propa¬ 
gated  between  measurement  epochs  as  follows: 


where 


(3.31) 


8t  =  ti  +  i-ti.  (3.32) 

The  two  frequencies  are  determined  by  one  of  the  previously  described  fre¬ 
quency  estimation  techniques. 

The  measurement  model  which  relates  the  displacement  vectors  to 
the  states  is  inherently  nonlinear  as  the  states  are  embedded  in  the  ele¬ 
ments  of  the  attitude  matrix.  The  nonlinear  measurement  equation  is 
given  by 


For  the  gain  and  covariance  update  calculations  the  linearized  measure¬ 
ment  connection  or  sensitivity  matrix  is 
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_  5h 
"dX 


(3.34) 


Simulated  data  is  generated  using  a  nutation  angle  of  5  deg,  an  iner¬ 
tia  ratio  of  1.3,  and  a  sample  frequency  of  2  Hz.  The  two  driving  frequen¬ 
cies,  cop  and  CO/,  are  -0.1561  rad/s  and  0.6789  rad/s,  respectively. 


3.5.3  Ground  Experiment  Batch  Processor 

The  batch  processing  algorithm  for  the  experimental  data  is  very 
similar  to  the  one  employed  for  the  computer  simulation.  The  estimated 
values  are  the  angular  position  of  the  antenna  about  the  body  zB-axis,  y, 
the  amplitude,  0,  of  the  tipping  motion,  the  argument  for  the  equation  for 
a(t),  X,  and  the  two  angles,  a  and  rj.  The  resulting  state  vector  is 

X  =  [y  0  X  rj  a]  •  (3.35) 

As  described  in  section  3.3.2,  the  induced  motion  is  designed  to  sim¬ 
ulate  torque  free  satellite  motion.  As  a  result,  the  rj  and  a  stay  fixed.  The 
magnitude  and  rate  of  the  tipping  motion  remain  constant.  The  position 
angles,  y  and  X,  change  according  to  the  two  driving  rates,  co  and  coa, 
respectively.  The  transition  between  measurement  times  is  therefore 
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co8t 


X(ti  +  1)  =  X(tj)  + 


0 

®a8t 


0 

0 


(3.36) 


The  driving  frequencies  are  estimated  using  the  frequency  analysis  tech¬ 
niques.  The  observations  equations  are  the  same  as  shown  in  Equations 
(3.33)  and  (3.34). 

The  phase  difference  A<|)  and  LOS  measurements  are  collected  at  a 
sample  rate  of  2  Hz  using  a  Trimble  TANS  Vector  receiver.  The  spin  and 
tipping  rates,  go  and  coa,  are  set  at  0.535  rad/s  and  0.597  rad/s,  respectively. 
The  tipping  amplitude  is  5.3  deg. 


3.5.4  Estimation  Results 

The  estimation  errors  for  rj  and  a  are  approximately  the  same  when 
using  either  the  cross  product  or  batch  algorithms.  When  N  is  400  or 
greater,  the  errors  are  under  0.1  deg.  The  batch  processor  allows  for  the 
estimation  of  additional  parameters  as  well.  For  the  largest  batch  sizes 
the  simulation  nutation  angle  and  the  experimental  tip  angle  are  both 
determined  to  within  0.1  deg.  Estimation  of  the  phase  angles  is  the  least 
accurate.  The  simulation  angles,  \| /  and  cp,  are  determined  to  about  0.2  deg. 
However,  the  experimental  angle,  y,  is  only  estimated  to  within  2  deg. 
Tables  3.3  and  3.4  shows  the  performance  of  the  attitude  methods  for  vari¬ 
ous  values  of  N  for  both  the  simulated  and  experimental  cases.  All  of  the 
data  is  taken  at  2  Hz.  For  the  simulation,  the  frequency  is  5  RPM,  and  the 
nutation  angle  is  5  deg.  For  the  experiment,  the  spin  frequency  is  5  RPM 
and  the  tipping  rate  is  2  deg/s. 
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Table  3.3.  Performance  of  attitude  methods  on  simulated  data. 


Method 

Angle 

True  value 

Estimated  value  (deg) 

(deg) 

N  =  200 

N  =  400 

N  =  800 

dr  Cross- 

T1 

63.43 

63.46 

63.46 

63.46 

Product 

o 

65.91 

65.90 

65.89 

65.91 

*n 

63.43 

63.47 

63.45 

63.45 

a 

65.91 

65.93 

65.92 

65.91 

Batch  Processor 

e 

5.00 

4.96 

4.98 

4.99 

¥ 

90.00 

91.82 

90.34 

89.83 

<p 

90.00 

88.23 

89.59 

90.18 

Table  3.4.  Performance  of  attitude  methods  on  experimental  data. 


Method 

Angle 

Reference 

value 

Estimated  value  (deg) 

(deg) 

N  =  200 

N  =  400 

N  =  800 

dr  Cross- 

T| 

0.60 

0.34 

0.52 

0.56 

Product 

c 

0.07 

0.11 

0.07 

0.08 

n 

0.60 

0.51 

0.52 

0.59 

a 

0.07 

0.10 

0.07 

0.10 

Batch  Processor 

e 

5.32 

5.13 

5.19 

5.29 

Y 

-47.38 

-45.15 

-44.74 

-49.01 

a 

1.01 

0.78 

0.82 

1.11 

3.6  Conclusions 

Several  new  algorithms  for  using  GPS  measurements  to  determine 
the  attitude  and  angular  rates  of  a  spinning  satellite  are  presented  and 
demonstrated  using  both  simulated  orbital  data  and  ground  test  data.  The 
AR  method  is  the  best  algorithm  for  determining  spin  rates.  The  dr  cross- 
product  method  is  the  computationally  simplest  and  fastest  approach  for 
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estimating  the  angular  momentum  axis.  The  batch  processor,  however, 
produces  estimates  of  more  vehicle  parameters.  For  batch  sizes  of  200 
samples,  corresponding  to  100  seconds  of  data  for  a  vehicle  spinning  at  5 
RPM,  spin  rates  are  determined  to  within  0.01  rpm,  spin  axis  orientation 
is  determined  to  within  0.25  deg,  and  phase  to  within  2  deg.  Comparable 
performance  could  be  expected  for  spacecraft  in  low  earth  orbit,  thus  satis¬ 
fying  the  requirements  of  many  small  satellite  missions. 


CHAPTER 4 


SNR  AND  OFF-BORESIGHT  ANGLE  RELATIONSHIP 

4.1  Introduction 

This  chapter  explores  the  relationship  between  the  signal-to-noise 
ratio  (SNR)  of  a  GPS  measurement  and  the  orientation  of  the  antenna 
boresight  vector.  It  begins  by  defining  the  SNR  observable.  This  is  fol¬ 
lowed  by  an  analysis  of  several  factors  that  affect  this  measurement.  The 
most  important  is  the  gain  pattern  of  the  GPS  receiver  antenna.  Other 
influences  are  the  path  loss,  the  transmitting  antenna  gain  pattern,  the 
transmission  power,  and  multipath.  Two  methods  are  developed  for  quan¬ 
tifying  the  relationship  between  the  SNR  and  the  antenna  boresight.  The 
first  technique  uses  a  statistical  analysis  of  actual  data  to  create  a  map¬ 
ping  function  between  the  SNR  measurements  and  the  angles  between  the 
boresight  and  the  LOS  vectors.  These  angles  are  called  the  off-boresight 
angles  and  are  denoted  by  the  symbol  a.  The  second  creates  a  mapping 
function  from  an  antenna  gain  pattern  scaled  to  match  the  GPS  receiver 
SNR  values.  The  chapter  concludes  with  a  method  for  calibrating  the  ini¬ 
tial  alignment  of  the  boresight  vector  based  on  SNR  measurements  and 
attitude  information. 
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4.2  SNR  Measurements 

In  addition  to  code  and  carrier  phase  measurements,  all  GPS 
receivers  make  some  type  of  signal-to-noise  ratio  (SNR)  measurement;  pri¬ 
marily  as  a  basis  for  checking  signal  integrity  and  maintaining  signal 
tracking.  Measurements  with  a  very  low  SNR  are  considered  unreliable 
and  may  not  be  used  in  a  navigation  solution.  Many  receivers  provide  the 
SNR  or  power  measurement  in  terms  of  the  carrier-to-noise  spectral  den¬ 
sity  (C/N0)  measured  in  dB-Hz.  In  the  Trimble  TANS  receivers  it  is  repre¬ 
sented  in  amplitude  measurement  units  (AMU).  The  relationship  between 
C/N0  and  AMU’s  is  given  by  [Cohen,  1992] 

C/N0  (dB-Hz)  =  101og10(BW*AMU2)  (4.1) 

4.3  Receiver  Antenna  Gain  Pattern 

The  relationship  between  the  SNR  of  a  measurement  and  the 
antenna  boresight  depends  on  the  angle,  a,  between  the  boresight  vector 
and  the  line-of-sight  to  the  GPS  satellite,  and  the  gain  pattern  of  the  GPS 
receiver  antenna.  Figures  4.1  through  4.4  show  the  gain  patterns  of  four 
different  GPS  antennas  used  in  this  research.  The  plots  represent  the  gain 
variations  with  respect  to  a  for  a  single  azimuth  angle  defined  as  0  deg.  As 
demonstrated  in  these  figures,  typical  GPS  antennas  have  the  highest  gain 
along  the  boresight  vector  with  decreasing  gain  down  to  90  deg  off-bore- 
sight.  One  exception,  as  shown  in  Figure  4.2,  is  for  the  TANS  antenna 
which  has  a  dip  near  0  deg.  As  illustrated  in  Figure  4.5,  the  effect  of  the 
gain  pattern  causes  the  SNR  to  be  highest  for  small  a  angles  and  lowest  for 
large  angles.  The  azimuthal  variations  for  all  of  these  antennas  are  quite 
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Figure  4.4 
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small.  For  example,  the  TANS  antenna  varies  less  than  1  dB  at  small  a’s 
and  less  than  3  dB  at  a’s  near  90. 


Antenna  Boresight 


Figure  4.5.  SNR  and  Boresight  Relationship. 

A  method  to  estimate  the  boresight  attitude,  described  in  Chapter  5, 
is  based  upon  the  SNR  measurement  being  dependent  solely  upon  the 
value  of  a.  However,  several  other  factors  affect  the  SNR  such  as  the  GPS 
satellite  transmit  power,  the  antenna  pattern  of  the  GPS  transmitter, 
space  loss,  atmospheric  attenuation,  multipath,  and  receiver  tracking 
error.  While  a  number  of  these  are  not  predictable  and  must  be  treated  as 
noise,  adjustments  can  be  made  for  two  of  the  factors. 

4.4  Space  Loss 

Space  loss  is  defined  as  the  reduction  in  signal  power  due  to  propa¬ 
gation  through  free  space.  It  is  computed  by 
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(4.2) 


where  X  is  the  wavelength  of  the  transmitted  signal  and  p  is  the  propaga¬ 
tion  distance.  If  a  signal  is  transmitted  with  power,  P,  then  the  received 
power,  Pn  is 


(4.3) 


In  dB  this  equation  is  written  as 

Pr[dBW]  =  101og10P[W]  -Ls(dB)  , 


(4.4) 


where 

Ls(dB)  =  201ogloP-201og10(^).  (4.5) 

The  effect  of  space  loss,  as  evident  in  Equation  (4.3),  is  that  the 
received  power  from  the  most  distant  GPS  satellite  is  reduced  more  than 
those  from  the  closer  satellites.  This  is  an  SNR  relationship  that  is  inde¬ 
pendent  of  the  antenna  orientation.  To  remove  this  effect,  the  SNR  values 
are  adjusted  so  excessive  space  loss  is  added  back  into  the  measurement. 
The  reference  space  loss  is  defined  as  the  amount  of  signal  power  reduction 
at  a  reference  distance,  pref.  For  this  research  pref  is  set  to  19,000  km.  This 
value  is  chosen  because  it  is  approximately  equal  to  the  minimum  distance 
between  a  low  earth  orbiting  (LEO)  spacecraft  at  an  altitude  of  1000  km 
and  a  GPS  satellite  at  zenith.  Using  Equation  (4.4)  the  power  received, 
Pref,  at  the  reference  distance  is 
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P«t[dBW]  =  101ogloP[W]  +201oglo(|;)-201og10p„f.  (4.6) 

The  space  loss  at  the  reference  distance  is  182  dB.  Subtracting  Equation 
(4.4)  from  (4.6)  yields 

Pref  -  Pr  =  201ogloPr-201°glOPref  *  (4  7) 

The  resulting  difference  in  Equation  (4.7)  is  the  adjustment  that  needs  to 
be  applied  to  each  SNR  measurement.  If  the  LEO  satellite  is  orbiting  at 
the  reference  altitude  or  below,  then  every  SNR  adjustment  increases  the 
actual  measurement.  For  satellites  above  the  reference  altitude,  the  SNR 
adjustment  decreases  the  measurement  for  small  zenith  angles  and 
increases  it  for  large  zenith  ones. 

The  effect  of  the  SNR  adjustment  is  most  evident  when  the  bore- 
sight  vector  is  pointing  more  towards  the  Earth’s  horizon.  As  illustrated  in 
Figure  4.6,  GPS  satellite  #2,  near  the  horizon,  is  farther  away  than  #1 
which  is  closer  to  zenith.  Satellite  #1,  however,  is  at  a  much  larger  a  with 
respect  to  the  antenna  boresight.  While  the  SNR  for  satellite  #2  should  be 
much  larger  than  from  satellite  #1  based  upon  the  angular  relationship, 
the  greater  level  of  space  loss  reduces  this  difference.  However,  knowing 
the  transmission  distances  allows  for  the  compensation  of  the  space  loss 
and  thus,  the  preservation  of  the  relationship  between  a  and  SNR.  When 
the  boresight  is  pointing  in  the  zenith  direction,  this  effect  is  not  as  appar¬ 
ent.  Measurements  from  GPS  satellites  at  zenith  have  the  highest  SNR’s 
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due  to  the  a  relationship  and  they  also  have  the  smallest  space  loss.  Con¬ 
versely,  signals  from  GPS  satellites  near  the  Earth’s  horizon  have  the 
smallest  SNR’s  and  they  have  the  greatest  space  loss.  As  a  result,  the 
space  loss  effect  actually  strengthens  the  a  dependency. 

Sat.  #1 


4.5  Transmitting  Power  Gain  Pattern  Correction 

The  gain  pattern  of  the  transmitting  antennas  on  the  GPS  satellites 
is  designed  to  provide  an  approximately  uniform  signal  power  reaching  all 
areas  on  the  surface  of  the  earth.  As  shown  in  Figure  4.7.,  the  antenna 
gain  is  highest  at  the  outer  edge  of  the  beam  width  to  compensate  for  the 
greatest  space  loss  on  the  ground  and  lower  in  the  middle  of  the  pattern 
where  the  space  loss  is  lowest.  The  gain  drops  sharply  beyond  15  deg  off- 
boresight,  where  there  are  no  users  on  the  surface.  For  space  based  users, 
however,  the  space  loss  compensation  requirements  are  different.  As  a 
result,  the  transmitting  gain  pattern  does  not  exactly  negate  the  effects  of 
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Figure  4.7.  GPS  satellite  antenna  gain  pattern.  [Czopek  and  Shollen- 
berger,  1993] 

space  loss.  As  illustrated  in  Figure  4.8,  the  values  of  a  between  the 
receiver  antenna  and  the  two  GPS  satellites  are  approximately  equal. 
However,  the  signal  from  satellite  #2  is  transmitted  at  a  larger  angle  in  the 
GPS  satellite  gain  pattern  as  compared  with  satellite  #1.  Thus,  it  is  an 
additional  source  for  the  variation  of  received  signal  power  that  is  indepen¬ 
dent  of  a. 

To  find  the  signal  variation  due  to  the  transmitting  antenna  gain 
pattern,  the  angle  between  the  LOS  vector,  e,  and  the  GPS  satellite 
antenna  boresight,  bops  >  must  be  computed  as  illustrated  in  Figure  4.8. 
The  boresight  of  the  GPS  satellite  antenna  is  approximately  in  the  direc¬ 
tion  of  the  Earth’s  center  and  can  be  computed  by 
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Sat.  #2 


Figure  4.8.  Angle  in  transmitting  antenna  pattern. 


I*  GPS  =  - 


RGPS 

rgps|  ’ 


(4.8) 


where  RGPS  is  the  position  of  the  GPS  satellite.  The  off-boresight  angle  of 
the  receiver  with  respect  to  the  GPS  satellite  boresight  is 


aGPS  =  cos1  (-e  •  bGps)  • 


(4.9) 


Figure  4.9  shows  an  example  of  the  power  level  adjustments  made 
for  a  LEO  satellite.  The  satellite  is  in  a  circular  equatorial  orbit  at  an  alti¬ 
tude  of  1000  km.  GPS  satellite  #32  (PRN  #1)  is  being  tracked.  The  plot 
covers  a  period  from  when  the  GPS  satellite  first  enters  the  field  of  view  of 
the  antenna  to  when  it  exits.  At  the  start  of  the  track,  the  LEO  and  GPS 
satellites  are  farthest  apart.  They  approach  each  other  until  some  mini¬ 
mum  distance  is  reached  and  then  grow  steadily  farther  apart.  The  thin 
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line  is  for  the  SNR  adjustment  for  the  space  loss.  As  is  expected,  the  com¬ 
pensation  is  the  highest  at  the  beginning  and  end  of  the  tracking.  The  dot 
and  dash  line  is  for  the  pattern  correction  adjustment.  When  the  GPS  sat¬ 
ellite  is  the  closest,  the  line-of-sight  is  near  the  center  of  the  pattern.  As  a 
result,  the  adjustment  is  close  to  zero.  At  two  points  in  the  track,  the  angle 
between  the  LOS  vector  and  the  GPS  satellite  antenna  boresight  vector 
reaches  about  10  deg.  These  points  correspond  to  the  greatest  gains  in  Fig¬ 
ure  4.7  and  therefore,  are  the  points  where  the  SNR  is  adjusted  downward 
the  most.  The  thick  line  is  the  combined  SNR  adjustment.  The  maximum 
increase  in  SNR  is  around  2.5  AMU  and  the  maximum  decrease  is  about  -2 
AMU.  These  values  are  about  5  percent  of  a  typical  total  SNR  measure¬ 
ment.  In  Figure  4.10  the  combined  adjustment  is  shown  as  a  function  of 
the  zenith  angle  between  the  LEO  satellite  and  the  GPS  satellite.  The 
greatest  increases  occur  between  65  and  90  deg.  The  largest  decreases  are 
between  20  and  60  degrees. 

4.6  Multipath  Effects 

Multipath  is  the  effect  caused  when  an  electromagnetical  signal 
reaches  an  antenna  by  an  indirect  path  through  one  or  more  reflections  as 
well  as  through  the  direct  path  [Comp,  1996].  Multipath  causes  the  SNR 
measurements  to  have  periodic  fluctuations.  The  frequencies  of  the  fluctu¬ 
ations  are  dependent  upon  the  distances  between  the  antenna  and  the 
reflectors  and  the  satellite  motion.  Figure  4.11  shows  an  example  of  how 
multipath  influences  the  SNR  measurement  over  time.  The  data  is  taken 
with  a  ground  based  antenna  from  a  single  GPS  satellite.  The  elevation 
angle  between  the  GPS  satellite  and  the  antenna  is  increasing  over  the 
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Figure  4.11.  The  effect  of  multipath  on  SNR. 

span  of  the  data.  The  SNR  increases  because  of  the  decrease  in  a.  How¬ 
ever,  there  are  also  rapid  fluctuations  due  to  multipath  from  distant 
objects.  Comp  [Comp,  1996]  develops  methods  for  estimating  multipath 
contributions  for  the  purpose  of  removing  phase  variations  due  to  its  effect 
[Axelrad  et  al.,  1996] .  The  multipath  effect  tends  to  be  a  strong  function  of 
the  LOS.  While  it  is  possible  to  construct  a  correction  map  based  on  exten¬ 
sive  data  sets,  the  technique  is  computationally  intensive  and  is  not  neces¬ 
sary  for  the  level  of  accuracy  achievable  in  this  work.  As  a  result, 
multipath  is  not  an  effect  on  SNR  that  can  be  removed  in  the  same  manner 
as  for  the  space  loss  and  pattern  correction.  Instead,  it  is  treated  as  mea 


surement  noise. 
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4.7  Variations  Due  to  Individual  GPS  Satellites 

The  GPS  satellites  are  designed  so  that  the  received  signal  power  on 
the  surface  of  the  Earth  does  not  exceed  -153.0  dBw.  The  signal  power  is 
generally  highest  when  the  satellite  is  new  and  drops  off  as  the  satellite 
approaches  its  end  of  life.  Over  the  lifetime  of  a  satellite,  the  power  loss  is 
expected  to  be  less  than  7.0  dB  [Ward,  P.,  1996].  Several  variables  affect 
the  power  output  when  compared  among  satellites.  Each  satellite  is  in  a 
different  stage  of  their  design  lifetime.  Each  satellite  receives  a  different 
level  of  solar  flux  and  experiences  different  environmental  conditions. 
There  are  also  variations  among  the  physical  characteristics  of  each  satel¬ 
lite.  The  most  evident  physical  difference  is  caused  by  having  three  differ¬ 
ent  satellite  models  currently  in  operation. 

These  power  variations  do  not  alter  the  a  versus  SNR  dependence 
like  the  space  loss  and  transmitting  gain  pattern  effects.  They  are  inde¬ 
pendent  of  the  receiver  type,  position,  and  orientation.  Their  influence  is 
apparent  on  a  satellite  by  satellite  basis.  For  these  variations  to  be  signif¬ 
icant,  a  plot  of  a  versus  SNR  for  an  individual  satellite  needs  to  show  a 
bias  when  compared  to  a  plot  of  a  versus  SNR  for  the  entire  constellation. 
In  other  words,  over  the  entire  range  of  SNR  values,  the  average  a  mea¬ 
surement  for  each  SNR  is  either  consistently  higher  or  consistently  lower 
than  the  average  for  the  entire  constellation. 

Data  from  the  CRISTA-SPAS  satellite,  discussed  more  thoroughly  in 
Chapter  6,  does  not  demonstrate  a  significant  bias.  This  set  is  chosen 
because  there  is  sufficient  data  covering  almost  every  GPS  satellite.  Fig¬ 
ures  4.12  through  4.15  show  four  examples  of  a  versus  SNR  covering  both 
block  II  and  IIA  GPS  satellites  and  a  range  of  satellite  lifetimes.  While  it 
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is  evident  that  the  average  a  values  from  the  individual  satellite  are  differ¬ 
ent  than  the  constellation  average,  the  differences  are  not  biased  towards 
either  higher  or  lower  values.  The  individual  averages  fluctuate  about  the 
constellation  average,  remaining  within  the  lcr  bounds. 

4.8  SNR  to  a  Mapping  Function 

The  key  to  a  reliable  attitude  solution  is  a  valid  mapping  of  SNR  to 
a.  As  is  discussed  in  Chapter  6  for  some  actual  GPS  data,  the  true  rela¬ 
tionship  between  SNR  and  a,  is  far  from  one-to-one.  The  ideal  relationship 
presented  in  section  4.3  does  not  accurately  reflect  the  variations  shown  in 
real  data.  An  adequate  mapping  function,  however,  can  still  be  estab¬ 
lished.  In  this  section  two  basic  approaches  to  create  this  function  are 
investigated-  the  first  creates  a  mapping  based  on  a  calibration  data  set 
and  the  second  uses  measured  antenna  gain  patterns.  The  latter  is  similar 
to  suggestions  in  [Hashida  and  Unwin,  1993]  and  [Serrano,  et  al.,  1995]. 

4.8.1  Calibration  Based  Mapping  Function 

A  calibration  data  set  consists  of  a  large  set  of  LOS  vectors,  mea¬ 
sured  SNR’s,  and  accurately  measured  a’s.  A  sufficiently  large  set  is  one 
that  includes  the  full  range  of  possible  viewing  angles  to  each  GPS  satel¬ 
lite.  It  requires  the  antenna  to  be  in  a  known  orientation  for  accurate  com¬ 
putation  of  a.  It  should  be  stressed  that  since  the  SNR  to  a  relationship  is 
not  one-to-one,  each  SNR  has  a  range  of  corresponding  a  values  associated 
with  it. 

Given  an  appropriate  calibration  data  set,  the  mapping  function  is 
generated  as  follows: 

1)  divide  the  SNR  values  into  discrete  bins; 
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2)  collect  the  accurate  a’s  corresponding  to  the  SNR  values  in 
each  bin; 

3)  compute  the  mean  a  and  standard  deviation  for  each  bin. 

In  applying  the  mapping  to  subsequent  data,  the  mean  a  is  the  best  esti¬ 
mate  for  the  bin.  In  other  words,  it  is  the  value  of  a  for  all  SNR’s  in  the 
bin.  The  standard  deviation  of  the  a’s  determines  the  appropriate  weight¬ 
ing.  To  quantify  the  overall  quality  of  the  mapping  function  the  value  Ap  is 
computed.  It  is  defined  as  the  mean  of  the  residuals  between  the  data 
points  and  a  5th  order  polynomial  fit  to  the  points. 

Calibration  data  sets  may  be  established  on  the  ground  or  in  orbit. 
In  this  research,  calibration  from  the  actual  data  being  processed  for  esti¬ 
mation  is  used  to  quantify  the  best  possible  performance.  This  is  called  the 
reference  result. 

The  best  practical  calibration  source  is  from  actual  orbit  data  from 
the  on-board  GPS  receiver.  During  times  when  the  attitude  is  known  from 
an  external  reference,  the  a’s  are  computed  and  SNR  data  are  recorded. 
Using  this  information,  a  mapping  function  is  created.  This  method  has 
the  advantage  that  the  mapping  function  is  based  on  data  taken  under  the 
same  conditions  as  the  present  data.  The  same  antenna  is  used,  the  noise 
in  the  receiver  is  similar,  the  line  biases  are  identical,  and  the  multipath 
environment  is  the  same.  Ideally,  the  calibration  data  include  measure¬ 
ments  taken  from  all  visible  GPS  satellites  throughout  the  orbits.  The 
major  disadvantage  is  that  SNR-based  attitude  estimation  is  not  available 
immediately.  One  scenario  where  this  might  be  useful  is  if  the  nominal 
attitude  reference  sensor  fails,  a  GPS  receiver  used  for  positioning  infor¬ 
mation  could  then  be  employed  for  attitude  determination.  It  might  also 
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be  used  for  faster  integer  ambiguity  resolution  in  a  GPS  phase-based  atti¬ 
tude  algorithm. 

The  second  calibration  source  is  from  data  collected  prior  to  the 
launch  of  the  satellite.  This  is  obtained  from  either  ground  tests  or  from  a 
previous  satellite  mission.  A  major  advantage  of  this  method  is  that  the 
mapping  function  is  available  for  immediate  estimation.  Unlike  on  orbit 
calibration,  however,  several  factors  could  be  significantly  different 
between  the  calibration  data  and  the  actual  measurements.  If  data  is  used 
from  a  previous  satellite  mission,  the  multipath  and  noise  environment 
could  be  significantly  different.  Also,  the  visibility  conditions  might  not  be 
the  same.  On  the  ground,  data  collection  would  most  likely  occur  on  a 
spacecraft  mock-up.  Multipath  from  the  ground  or  other  surroundings 
would  be  added  resulting  in  a  different  SNR  behavior.  The  ground  receiver 
would  see  a  different  set  of  GPS  satellites  and  it  would  be  moving  much 
slower  in  relation  to  them.  Finally,  the  environmental  conditions  of  actual 
receiver  operation  would  be  different. 

The  SNR  to  a  mapping  can  be  readily  updated  whenever  SNR  data 
and  an  external  reference  are  available.  This  could  be  important  for  a  sat¬ 
ellite  that  uses  a  number  of  distinct  orientations  over  its  mission  lifetime. 

4.8.2  Gain  Pattern  Based  Mapping  Function 

An  alternative  method  for  generating  a  mapping  function,  which 
does  not  depend  upon  a  calibration  data  set,  involves  using  a  gain  pattern 
of  the  antenna.  Such  patterns  are  routinely  measured  in  an  anechoic 
chamber.  The  gain  pattern  values  must  be  scaled  by  the  nominal  received 
power  in  order  to  match  the  actual  output  values  of  the  GPS  receiver.  For 
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example,  the  top  graph  in  Figure  4.16  shows  the  measured  gain  pattern  in 
dB  of  a  patch  antenna.  The  bottom  graph  shows  the  measured  gain  pattern 
scaled  to  match  the  SNR  in  AMU’S  of  a  TANS  Vector  receiver.  The  conver¬ 
sion  equation  is 

AMU  =  ^antenna  7^0’  <410> 

where  P0  is  the  arriving  signal  strength. 


Gain  (db) 


Figure  4. 16.  Scaling  of  measured  gain  pattern  values  to  receiver 
output  values. 


The  scaling  factors  may  be  derived  by  monitoring  the  SNR  output  of 
the  receiver.  It  is  assumed  that  over  a  20  minute  period  one  or  more  of  the 
GPS  satellites  has  a  maximum  a  of  80  to  90  deg.  These  should  correspond 
to  the  lowest  measured  SNR’s.  Similarly,  one  or  more  GPS  satellites 
should  be  within  5  deg  of  the  boresight  corresponding  to  the  highest  SNR’s. 
For  example,  in  the  GANE  data  set,  discussed  thoroughly  in  Chapter  6,  the 
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highest  a  in  a  typical  20  minutes  segment  of  data  is  80  deg.  This  measure¬ 
ment  has  a  SNR  of  4.5  AMU.  For  the  entire  data  set,  about  1%  of  the  SNR 
values  are  below  4.5.  The  smallest  a  is  4  deg.  This  measurement  has  an 
SNR  of  29  AMU.  About  1%  of  the  SNR  values  are  above  this  number. 
Assuming  that  these  extreme  values  from  the  sample  20  minutes  of  data 
are  good  representatives  of  any  block  of  data,  they  are  used  to  create  scal¬ 
ing  factors  for  the  gain  pattern. 

The  scaling  factors,  k,  are  computed  from  the  ratio 


where  Kr  is  the  SNR  in  AMU’S  from  the  receiver  at  one  of  the  extreme 
angles  and  Kg  is  the  corresponding  pattern  gain  in  dB  at  this  angle.  Kr  is 
computed  from  the  mean  of  the  expected  range  of  SNR’s  for  this  angle.  For 
example,  in  the  GANE  data  the  la  variation  in  the  SNR’s  corresponding  to 
80  deg  cover  the  lowest  2%  of  the  recorded  SNR  measurements .  As  a 
result,  Kr  is  computed  from  the  mean  of  this  range.  As  shown  in  Figure 
4.16,  the  lower  gain  pattern  values,  plotted  in  the  top  graph,  are  scaled  by 
the  value  of  k  for  the  lowest  SNR’s,  while  the  higher  gain  pattern  values 
are  scaled  by  k  for  the  highest  SNR’s.  The  mid  range  values  are  scaled  by 
an  average  of  these  two  factors. 

The  greatest  advantage  to  this  method  is  that  the  mapping  function 
is  easy  to  generate  and  is  available  for  immediate  use  on  orbit.  However, 
the  environment  under  which  the  SNR  mapping  is  generated  is  the  most 
different,  does  not  include  multipath,  and  can  not  account  for  GPS  satellite 
variations.  It  also  may  not  work  well  for  satellites  in  certain  orientations 
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for  which  the  entire  a  range  is  not  covered  in  a  short  time  period.  The  scal¬ 
ing  factors,  in  this  case,  can  be  easily  adjusted  as  more  data  becomes  avail¬ 
able. 

4.8.3  Mapping  Function  Summary 

Of  all  the  methods  to  generate  a  mapping  function,  on-orbit  calibra¬ 
tion  is  the  best  choice  in  terms  of  accuracy.  The  conditions  under  which  the 
function  is  computed  remain  consistent  throughout  the  mission.  There¬ 
fore,  the  mapping  function  data  will  correspond  very  closely  to  the  behav¬ 
ior  of  the  measurements  used  for  the  actual  axis  estimation. 

After  generating  this  function  by  either  the  data  calibration  or  gain 
pattern  methods,  attitude  estimation  can  occur.  During  each  estimation 
epoch,  one  SNR  value  is  measured  for  each  GPS  satellite  in  view.  Each  of 
the  corresponding  bin  numbers  are  found.  In  those  bins  the  associated  a 
and  weight  are  located.  Chapter  5  describes  the  process  by  which  the 
angles,  measurement  weights,  and  LOS  vectors  are  used.  The  result  is  an 
instantaneous  estimate  of  the  antenna  boresight  direction. 

4.9  Estimation  of  the  Boresight  Alignment 

To  further  refine  the  mapping  function,  an  auxiliary  calibration  step 
is  required  to  determine  the  offset  between  the  nominal  body  orientation  of 
the  boresight  vector  and  its  true  orientation.  The  nominal  boresight  vector 
in  the  body  frame  can  be  defined  from  a  priori  knowledge  of  how  the 
antenna  is  mounted  on  the  spacecraft  structure.  For  example,  suppose  the 
spacecraft  body  frame  is  defined  so  that  one  face  is  called  the  xByB-plane.  If 
a  GPS  patch  antenna  is  mounted  with  its  plane  parallel  to  the  xByB -plane, 
then  the  boresight  vector  is  along  the  zB-axis.  If  the  patch  antenna  is  not 
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perfectly  parallel  to  this  plane,  then  there  is  some  constant  misalignment 
as  shown  in  Figure  4.17.  Another  possibility  is  if  the  gain  pattern  of  the 
antenna  is  not  maximum  along  the  nominal  boresight  axis.  The  offset 
must  be  determined  from  a  collection  of  GPS  data  and  a  knowledge  of  the 
true  attitude. 


As  discussed  in  section  4.8.1,  construction  of  a  calibration  data 
based  mapping  function  requires  knowledge  of  the  true  orientation  of  the 
antenna  boresight  vector  and  the  LOS  vectors.  By  knowing  the  true  atti¬ 
tude  of  the  spacecraft  with  respect  to  some  local  system,  the  LOS  vectors  in 
the  spacecraft  frame  are  computed  by  using  the  rotation 

ef  =  BCL§f.  (4.12) 

Each  oq  in  the  calibration  set  is  then  determined  from  the  relationship 
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(6B)  ftB  =  cosaj.  (4.13) 

The  batch  of  oq  values  is  used  to  estimate  the  calibration  function,  F,  such 
that 

F  (SNRj)  =  cosaj .  (4.14) 

B 

However,  if  the  true  orientation  of  the  antenna  in  the  body  frame,  n  ,  is 
not  available  due  to  some  unknown  misalignment,  then  the  true  values  of 
a  cannot  be  determined  and  thus,  F,  cannot  be  estimated. 

While  F  is  unknown,  other  information  is  available  based  on  the 
assumption  that  Equation  (4. 14)  is  valid.  Every  LOS  vector  satisfying  the 
conditions  that 

e^e-  (4.15) 

and 

SNR^SNRj,  (4.16) 

implies  that 

ctj-cij  (4.17) 

This  leads  to  the  relationship  that 

F(SNRi)  -F(SNRj)  .  (4.18) 

Combining  Equations  (4.13),  (4.14),  and  (4.18),  an  equation  for  n  ,  involv¬ 
ing  the  difference  of  two  LOS  vectors  associated  with  approximately  equal 
SNR  measurements  is 
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(4.19) 


For  any  pair  of  LOS  vectors,  there  are  an  infinite  number  of  possible 
solutions  for  n  .  However,  the  actual  solution  must  simultaneously  sat¬ 
isfy  Equation  (4. 19)  for  all  pairs  of  LOS  vectors  meeting  the  conditions  set 

B 

in  Equations  (4.14)  and  (4.15).  This  leads  to  a  procedure  to  estimate  n  . 
The  steps  are  as  follows: 

1)  group  all  of  the  SNR  measurements  into  discrete  bins; 

2)  group  the  LOS  vectors  according  to  the  bins  of  their  corre¬ 
sponding  SNR  measurements; 

3)  for  each  bin,  k,  form  the  matrix 


B  T 

(e2)  - 

,.B  T 
(ej) 

Gk  = 

(®3)T  - 

B  T 

(eB) 

(4)T- 

B  T 

(4) 

(4.20) 


where  N  is  the  number  of  measurements  in  the  bin; 

4)  combine  Equations  (4.19)  and  (4.20)  to  yield 


Gi 

G2 


=  0nB~O, 


where  M  is  the  number  of  SNR  bins; 

5)  compute  nB . 


(4.21) 
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In  the  ideal  case,  n  spans  the  null  space  of  g,  however,  for  real 

data  it  is  the  vector  that  produces  the  smallest  norm  when  g  is  multiplied 

B 

by  it  [Press,  et  al.,  1992].  n  can  be  computed  by  first  performing  a  singu¬ 
lar  value  decomposition  on  g  to  express  this  matrix  as  a  product  of  a  col¬ 
umn-orthogonal  matrix,  U,  a  diagonal  matrix,  S,  and  a  orthogonal  matrix, 
V,  such  that 

g  =  USVT.  (4.22) 

The  diagonals  of  S  are  the  singular  values  and  are  equal  to  the  norms  of 

the  columns  of  g\ .  If  g  is  truly  singular,  then  one  of  the  singular  values 

B 

is  exactly  zero  and  n  is  the  corresponding  column  of  V.  Because  of  the 

approximations  involved  in  generating  g,  there  is  no  singular  value  equal 

B 

to  zero.  Therefore,  n  is  approximated  by  the  column  of  V  associated  with 

B 

the  smallest  one.  To  reliably  estimate  n  a  sufficiently  large  set  of  SNR 
and  LOS  measurements  is  required  as  well  as  an  accurate  estimate  of  the 
vehicle  attitude  over  the  duration  of  the  data.  For  this  research,  the  entire 
span  of  data  is  used  for  each  set.  Further  work  is  necessary  to  determine 
the  set  size  requirements  for  a  given  level  of  accuracy. 

It  is  interesting  to  note  that  in  the  case  of  a  constant  attitude  vehi¬ 
cle,  the  solution  to  Equation  (4.21)  is  an  estimation  of  the  boresight  vector 
in  the  local  frame.  By  defining  the  body  frame  to  be  aligned  with  the  local 
frame,  the  boresight  vector  alignment  in  the  body  frame  is  also  its  orienta¬ 
tion  with  respect  to  the  known  local  system.  However,  multiple  LOS  vec¬ 
tors  for  each  SNR  bin  are  required  for  Equation  (4.21)  to  be  solvable.  This 
is  highly  unlikely  for  any  single  measurement  epoch.  Therefore,  Equation 
(4.21)  can  only  be  solved  over  a  batch  of  measurement  times.  When  the 
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vehicle  attitude  is  not  constant,  the  solution  for  Equation  (4.21)  can  only  be 
used  to  estimate  the  boresight  vector  in  the  local  system  at  any  single  time. 
As  previously  stated,  one  measurement  tune  is  not  sufficient  for  this  com¬ 
putation.  Over  a  batch  of  measurement  epochs,  the  attitude  is  changing 
and  as  a  result,  the  alignment  of  the  vehicle  body  with  respect  to  the  local 
system  is  changing.  However,  the  alignment  between  the  boresight  vector 
in  the  body  system  and  any  other  axis  in  the  body  frame  is  constant. 
Therefore,  as  shown  in  the  procedure  outlined  in  the  last  paragraph,  the 
solution  for  Equation  (4.21)  can  always  be  used  to  estimate  this  offset. 

If  truth  attitude  is  not  available  and  the  mapping  function  is  created 
from  the  antenna  gain  pattern,  then  a  mechanical  drawing  of  the  antenna 
mounted  on  the  spacecraft  is  necessary.  An  estimate  of  the  boresight  vec¬ 
tor  is  possible  without  the  drawing.  However,  the  relationship  between  the 
antenna  orientation  and  the  orientation  of  any  instrument  on  the  vehicle 
remains  unknown. 

4.10  Conclusions 

The  relationship  between  the  SNR  measurement  from  a  GPS 
receiver  and  the  orientation  of  the  antenna  boresight  vector  is  based  pri¬ 
marily  on  the  gain  pattern  of  the  antenna.  LOS  vectors  to  GPS  satellites 
that  are  at  small  angles  with  respect  to  the  antenna  boresight  correspond 
to  the  highest  SNR  measurements.  LOS  vectors  that  make  large  angles 
with  the  boresight  vector  correspond  to  the  lowest  SNR  measurements. 
Other  factors,  such  as  space  loss  and  the  transmitting  antenna  gain  pat¬ 
tern  contribute  to  this  relationship  in  a  manner  that  must  be  compensated. 
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Power  fluctuations  from  individual  GPS  satellites  are  not  significant  from 
the  current  constellation.  Multipath  is  treated  as  measurement  noise. 

Because  this  relationship  for  actual  data  is  unknown  and  not  one-to- 
one,  a  mapping  function  is  created.  The  technique  used  in  this  research  is 
based  upon  using  either  a  calibration  data  set  or  an  antenna  gain  pattern. 
The  calibration  data  set  requires  SNR  data  from  an  antenna  in  a  known 
orientation.  The  gain  pattern  approach  requires  just  SNR  data.  Before 
creating  this  mapping  function,  the  true  boresight  alignment  in  the  body 
system  has  to  be  computed  through  an  initial  calibration  step.  After  this  is 
accomplished,  the  estimation  of  the  boresight  vector  can  proceed. 

In  the  next  chapter,  a  method  is  presented  that  uses  the  mapping 
functions  to  estimate  the  boresight  vector  of  an  antenna.  GPS  data  from 
several  sources  are  examined  in  Chapter  6.  The  data  are  processed  to  gen¬ 
erate  mapping  functions.  Chapter  7  presents  the  actual  estimation 
results. 


CHAPTER  5 


MAXIMUM  LIKELIHOOD  BASED  ESTIMATION  ALGORITHMS 

5.1  Introduction 

The  previous  chapter  details  the  relationship  between  the  SNR  mea¬ 
surements  and  a,  defined  to  be  the  angle  between  the  GPS  LOS  vectors 
and  the  orientation  of  the  antenna  boresight  vector.  Procedures  are  estab¬ 
lished  to  create  a  mapping  function  between  them.  In  the  first  part  of  this 
chapter  a  method  is  formulated  that  uses  a  and  SNR  as  observables  in  a 
maximum  likelihood  estimate  (MLE)  of  the  boresight  vector.  First  the 
geometry  of  the  problem  is  detailed.  Next,  the  mathematical  basis  of  the 
technique  and  its  application  to  the  SNR  measurements  are  presented. 
This  is  followed  by  a  derivation  of  the  estimation  error  covariance. 

The  second  part  of  this  chapter  outlines  the  application  of  this  single 
axis  approach  to  specific  attitude  estimation  problems.  Knowledge  of  a  sin¬ 
gle  boresight  vector  is  adequate  when  only  one  axis  is  required,  such  as  for 
the  spinning  satellite  case  in  which  the  orientation  of  the  angular  momen¬ 
tum  axis  is  needed.  For  this  case,  the  necessary  computations  under  vari¬ 
ous  antenna  mounting  positions  are  investigated.  In  the  case  of  a  three- 
axis  stabilized  satellite  a  minimum  of  two  non-colinear  axis  are  required. 
Unfortunately,  there  is  no  actual  satellite  data  available  for  either  of  these 
two  cases.  Performance  is  only  estimated  for  the  single-axis  results. 
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5.2  Estimation  Geometry 

A  standard  method  for  estimating  the  orientation  of  a  single  space¬ 
craft  axis  in  a  given  reference  frame  involves  two  angle  measurements 
between  the  unknown  axis,  n ,  and  two  known  vectors,  dj  and  e2 ,  referred 
to  in  this  research  as  the  reference  vectors.  As  shown  in  Figure  5.1,  each 
angle  measurement  constrains  the  solution  to  a  circle  on  the  surface  of  a 
unit  sphere,  centered  about  the  known  reference  vector.  The  radius  of  the 
circle  is  a  function  of  the  angle  measurement.  The  intersections  of  the  cir¬ 
cles  determine  the  only  possible  values  for  n .  With  perfect  instruments, 
the  addition  of  a  third  measurement  removes  the  ambiguity.  Any  addi¬ 
tional  measurements  also  provide  intersections  at  the  same  solution  point. 
This  is  provided  that  the  additional  reference  vectors  are  not  parallel  to 
the  others. 


Surface  of  Unit  Sphere 


Figure  5.1.  Single  axis  estimation  from  two  angle  measurements. 

In  many  small  satellite  missions,  the  reference  vectors  and  angle 
measurements  are  taken  from  an  Earth  horizon  sensor  and  sun  sensors. 
The  Earth  lines  are  determined  from  the  position  of  the  spacecraft  in  an 
Earth  centered  frame.  Because  the  Earth  is  very  close  to  the  spacecraft, 
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light  reflected  off  the  planet  is  not  considered  to  be  from  a  point  source. 
Instead,  a  horizon  sensor  is  used.  The  field-of-view  (FOV)  for  this  sensor 
sweeps  through  space  to  detect  two  rapid  changes  in  light  intensity 
between  space  and  the  Earth,  one  change  occurring  when  the  FOV  enters 
the  atmosphere  and  the  other  when  it  leaves.  Knowing  the  scanning  speed 
of  the  instrument,  the  time  between  horizon  detections,  and  the  size  of  the 
Earth,  the  angle  between  the  Earth  line  and  the  plane  of  the  sensors  is 
computed.  The  sun  lines  are  accurately  determined  from  the  known  posi¬ 
tion  of  the  spacecraft  with  respect  to  the  center  of  the  Earth  and  from  the 
known  position  of  the  Earth  with  respect  to  the  sun.  Because  the  sun  is  far 
away  from  the  spacecraft,  the  solar  rays  reaching  the  sensor  are  consid¬ 
ered  to  be  parallel  to  the  sun  line.  The  sun  sensor  measures  the  angle  of 
these  incoming  light  rays  with  respect  to  the  sensor  plane. 

GPS  also  may  be  used  in  the  same  manner  by  mapping  SNR  mea¬ 
surements  to  the  angle  a  formed  between  the  antenna  boresight  and  each 
GPS  LOS  vector,  e .  In  this  case,  the  antenna  boresight  is  the  unknown 
avis  to  be  determined  and  the  LOS  vectors  provide  the  reference  vectors  in 
an  appropriate  reference  frame.  This  is  diagramed  in  Figure  5.2.  LOS  vec¬ 
tors  are  determined  based  upon  the  known  positions  of  the  GPS  satellites 
and  the  user  spacecraft  carrying  a  GPS  receiver  and  antenna.  Like  the 
sun  and  Earth  lines,  the  LOS  vectors  are  assumed  to  be  highly  accurate. 
The  position  computations,  however,  do  not  need  to  be  very  accurate;  in 
fact  LEO  position  errors  of  up  to  350  km  only  contribute  as  much  as  1  deg 
error  in  the  LOS  computation. 

Assuming  a  perfect  relationship  between  SNR  and  a,  a  single  SNR 
measurement  from  one  satellite  gives  the  angle  between  the  unknown 


78 


n  M 


boresight  vector  and  the  LOS  vector.  In  three  dimensional  space,  this  a 
cone  of  possible  boresight  vectors  centered  on  the  LOS  vector.  On  the  unit 
sphere,  unit  vectors  are  represented  by  their  x,  y,  and  z  components  plotted 
on  the  surface.  For  the  case  of  one  LOS  vector  and  one  SNR  measurement, 
the  LOS  vector  is  a  single  point  on  the  surface  of  the  sphere  while  the  pos¬ 
sible  boresight  vectors,  corresponding  to  the  value  a  from  the  SNR  map¬ 
ping  function,  he  in  a  circle  about  this  point.  The  radius  of  the  circle  is  a 
function  of  a.  Figure  5.3  illustrates  both  the  three  dimensional  view  and 
the  unit  sphere  projection. 

An  additional  SNR  measurement  from  a  second  satellite  narrows 
the  possible  locations  of  the  boresight  vector  to  two  positions.  As  previ- 
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3D  view 


LOS  vector 


Top  view  of  unit  sphere 
* 


Possible  boresight  vectors 


Radius 
from  a 


Figure  5.3.  Boresight  Information  from  a  Single  GPS  Satellite. 

ously  stated,  the  individual  a  values  mapped  from  each  of  the  SNR  mea¬ 
surements  restrict  the  boresight  location  to  circles  about  the  LOS 
vectors.  The  true  boresight  location  must  he  on  both  of  these  circles.  This 
condition  occurs  at  the  intersections  between  the  two  LOS  centered  circles 
or  in  the  3-D  representation,  between  the  two  LOS  centered  cones,  as 
shown  in  Figure  5.4. 

3D  View  Top  view  of  unit  sphere 


Figure  5.4.  Multiple  SNR  and  LOS  vector  geometry. 
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Because  of  the  large  uncertainties  involved  in  the  actual  mapping 
functions,  the  GPS  case  is  better  illustrated  using  rings  instead  of  circles. 
The  intersection  of  two  rings  limits  the  possible  boresight  locations  to  two 
large  areas.  One  advantage  of  GPS  is  that  there  are  typically  as  many  as 
five  or  more  GPS  satellites  in  view  at  a  time.  Each  GPS  satellite  provides 
another  measurement  and  thus,  another  ring  of  possible  locations  for  n  on 
the  unit  sphere.  As  shown  in  Figure  5.5,  the  addition  of  more  rings  begins 
to  reduce  the  large  uncertainty  regions.  With  three  or  more  measurements 
an  optimal  estimate  of  n  can  be  determines  based  upon  all  the  information 
available.  The  optimal  choice  is  found  using  a  maximum  likelihood  esti¬ 
mate  (MLE). 


Single  SNR  Multiple  SNR 


Area  of  possible 
boresight  locations 


Figure  5.5.  Geometry  of  GPS  measurements. 
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5.3  Maximum  Likelihood  Estimates 

For  given  a  set  of  observed  values  Xj,  X2,  xn  taken  from  a  set  of 
random  variables,  X1(  X2,  Xn  and  a  set  of  unknown  parameters 
p2>  •••>  Pm>  maximum  likelihood  estimates  for  the  parameters  are  the 
values  that  maximize  the  joint  probabihty  density  function  [Devore,  1990] 

ffxi,  x2, xn;  pi,  P2, ...,  pm).  (51) 

The  MLE  is  desirable  because  it  is  consistent,  asymptotically  unbiased, 
asymptotically  normally  distributed  and  asymptotically  efficient.  In  addi¬ 
tion,  it  is  the  minimum  variance  imbiased  linear  estimate  [Maybeck, 
1994],  For  the  GPS-based  single  axis  estimation  case,  the  observed  values 
are  the  cosines  of  a  and  the  parameter  is  the  boresight  vector.  Therefore, 
the  problem  is  to  find  the  unit  vector  that  maximizes 

frcosa!,  cosa2, ...,  cosan;  n).  (5.2) 

The  weighted  least  squares  solution  is  one  unbiased  and  minimum 
variance  estimator  that  is  a  potential  MLE  for  (5.2).  The  boresight  vector 
is  the  solution  that  minimizes  the  cost  function 


N 


J(n)  =  \  X  ak  cosak-  (ek)  n 


(5.3) 


k  =  1 


where  ak  are  the  measurement  weights.  Shuster  [Shuster,  1989], [Shuster, 
1994]  shows  that  the  minimization  of  (5.3)  is  the  MLE  for  n  if  the  proba¬ 
bility  density  function  for  the  measurements,  cosak,  are  axially  symmetric 
about  the  n  and  if  the  weights  are  chosen  correctly.  The  correct  weights 
are  determined  to  be 
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where  are  the  variances  of  the  density  functions. 


(5.4) 


5.4  Shuster’s  Method 

In  this  research  the  axially  symmetric  distribution  assumption  is 
accurate  for  the  data  used  as  will  be  proven  in  Chapter  6.  Therefore,  the 
cost  function  to  be  minimized  is  given  by 


J(A) 


N 


-t 

-  y  — 

k  =  l°k 


cosak-  (ek) 


(5.5) 


where  the  optimal  solution  for  n  must  also  satisfy  the  unit  vector  con¬ 
straint.  The  form  of  Equation  (5.5)  is  expressed  for  any  reference  frame. 
The  representation  of  the  reference  vector  establishes  the  coordinate  frame 
in  which  the  estimate  for  A  is  resolved.  For  the  GPS  case,  the  most  conve¬ 
nient  and  useful  solution  is  expressed  in  the  local  level  reference  frame. 


5.4.1  Estimation  Procedure 

The  first  step  in  the  estimation  of  n  is  to  compute  an  initial  approx¬ 
imation,  no,  that  minimizes  Equation  (5.5)  without  the  unit  vector  con¬ 
straint.  This  solution  is  simply  the  weighted  least  squares  result  given  by 

n0  =  (HtWH)_1HtWc,  (5.6) 


where 
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and 


HT  .  [6l  s2 ...  eN] , 


(5.7) 
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(5.8) 


cT  =  ^coso^  cosa2  ...  cosa^j 


(5.9) 


The  result  is  then  normalized  to  produce  the  first  unit  vector  estimate 


(5.10) 

Because  nQ  is  the  vector  that  minimizes  Equation  (5.5)  in  the  most  general 
case,  the  value  of  the  cost  function  with  additional  constraints  on  the  solu¬ 
tion  can  never  be  lower.  As  a  result  the  value  of  Equation  (5.5)  at  nQ  rep¬ 
resents  a  lower  bound  for  the  actual  unit  vector  minimization. 

To  determine  the  optimal  estimate,  the  solution  vector  is  first 
expressed  in  the  terms  of  two  spherical  angles,  as  illustrated  in  Figure  5.6, 
such  that, 
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al  = 


sin<|>cos0 
sin  (|)  sin  0 
cos<j> 


(5.11) 


Figure  5.6.  Boresight  angles. 


Minimizing  the  cost  function  in  terms  of  these  two  angles  is  equivalent  to 
minimizing  it  subject  to  the  normal  constraint,  because  all  values  of  0  and 
<t>  produce  a  unit  vector  in  Equation  (5.11).  The  choice  of  the  angles  in  Fig¬ 
ure  5.6  is  arbitrary.  The  angle  <(>  can  be  defined  between  the  boresight  vec¬ 
tor  and  any  one  of  the  primary  axes.  However,  when  the  boresight  is 
parallel  to  the  primary  axis  from  which  $  is  defined,  <|)  is  zero  and  0  is  non¬ 
unique.  As  is  evident  in  Equation  (5.11),  the  x  and  y  components  of  n  are 
zero  for  any  value  of  0  when  <|)  is  zero.  Therefore,  a  definition  of  $  should  be 
chosen  that  is  not  close  to  zero. 

A  gradient  search  is  performed  on  the  two  angles  to  find  the  true 
minimum.  Making  the  definitions  that  0j  and  are  the  estimates  for  0  and 
<t>  at  iteration  i  and , 
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n 


and  8q 


(5.12) 


the  cost  function  written  in  terms  of  these  variables  and  expanded  in  a 
Taylor  series  is 


JOli)  =J(ni_i)+  8ni  + 


&li  +  — 


(5.13) 


Keeping  only  the  terms  up  to  those  quadratic  in  8q  ,  yields  the  solution  for 
the  correction  term 


(5.14) 


Defining  the  derivative  of  A  with  respect  to  the  two  angles,  <J)  and  0, 
as  the  matrix 


Sfi  _  8n  8n 
8t|  t  80  8(J) 


(5.15) 


the  first  derivative  of  J  with  respect  to  q  is  found  to  be 

|i(n)  =  -(|4)T  [HTWc  -  (HTWH)  n] ,  (5.16) 

where  H,  W,  and  c  are  the  matrices  defined  in  Equations  (5.7)  through 
(5.9),  respectively.  The  second  derivative  of  J  is 
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(tl) 
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[HTWc-  (HTWH)A] 


(5.17) 


The  second  term  in  Equation  (5.17)  is  much  smaller  than  the  first  term; 
therefore,  the  second  derivative  is  approximated  as 
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6rj  5r|T 
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(HTWH) 


(5.18) 


The  initial  estimates  for  0  and  (j)  are  determined  from  the  evaluation 
of  n0  and  Equation  (5.11).  The  correction  term  for  these  values  is  then 
computed  from  Equation  (5.14).  Updated  estimates  are  computed  from  the 
equation 
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(5.19) 


5.4.2  The  Shuster  Algorithm  Applied  to  SNR  Measurements 

For  the  GPS  case,  the  angles,  ak,  are  the  off-boresight  angles  derived 
from  the  SNR  measurements,  the  reference  vectors,  ek ,  are  the  LOS  vec¬ 
tors,  and  the  vector  to  be  estimated,  n ,  is  the  antenna  boresight.  The  com¬ 
putation  of  the  measurement  weights  is  taken  from  the  calibration 
function  as  described  in  Section  4.8.1. 

To  perform  instantaneous  attitude  estimation,  the  observations  are 
taken  simultaneously;  thus,  the  number  of  angles  is  limited  by  the  track¬ 
ing  capabilities  of  the  GPS  receiver.  In  other  words,  the  value  of  N  in 
Equation  (5.5)  is  equal  to  the  number  of  satellites  tracked.  It  is  interesting 
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to  note  that  the  unconstrained  measurement  gradient,  computed  by  Equa¬ 
tion  (5.7),  for  this  type  of  solution  is  merely  the  GPS  LOS  vector.  This 
means  that  an  indication  of  the  contribution  of  satellite  geometry  to  the 
attitude  is  described  by  the  position  dilution  of  precision  (PDOP)  generally 
used  to  quantify  navigation  performance.  This  is  computed  by 

PDOP  =  Vtrace  ( (HTH) -1) .  (5.20) 

Simpler  methods  for  computing  boresight  direction  based  upon  the 
angles  measurements  and  LOS  vectors  may  also  be  considered  for  compu¬ 
tational  reasons.  Two  techniques  evaluated  for  comparison  to  Shuster’s 
optimal  solution  are  a  simple  average  of  LOS  vectors  for  all  tracked  satel¬ 
lites,  and  the  n0  computed  from  normalized  result  of  Equation  (5.6)  used 
to  initialize  the  search  in  the  MLE  solution. 


5.4.3  Estimation  Error  Covariance 

Given  the  true  boresight  axis  and  assuming  the  prediction  errors 
are  normally  distributed  with  zero  mean,  the  theoretical  error  covariance 
matrix  based  on  the  MLE  solution  can  be  predicted  for  the  two  angles  in 
the  estimation  procedure  [Shuster,  19831.  The  inverse  of  the  covariance 
matrix  is 
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(5.21) 


Equation  (5.21)  gives  a  convenient  way  to  evaluate  the  expected  perfor¬ 
mance  of  the  boresight  estimation  procedure  on  the  test  data.  If  the  true 
attitude  is  known,  then  can  be  evaluated  from  the  values  of  n ,  H, 

and  W.  When  expanded,  the  covariance  matrix  is 
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E(502)  E  (805<t>) 
E  (505<t>)  E(5<))2) 


(5.22) 


Therefore,  the  expected  values  of  50  and  5(J)  are  simply  the  square  roots  of 
the  diagonals.  The  errors  computed  from  the  actual  results  should  agree 
with  these  expected  errors  to  a  level  of  certainty  of  la. 


5.5  Spinning  Satellite  SNR*Based  Attitude  Determination 

There  are  three  distinct  ways  a  GPS  antenna  might  be  mounted  on 
a  spinning  satellite.  The  preferred  orientation  is  with  the  plane  of  the 
antenna  lying  flat  on  one  of  the  satellite’s  faces  so  that  the  boresight  vector 
is  parallel  to  the  angular  momentum  axis.  The  second  is  with  the  plane  of 
the  antenna  mounted  on  the  side  of  the  satellite  so  that  the  boresight  vec¬ 
tor  is  perpendicular  to  the  angular  momentum  axis.  The  other  orientation 
is  with  the  boresight  mounted  between  these  two  extremes  at  an  acute 
angle  to  the  angular  momentum  axis. 


5.5.1  Boresight  Mounted  Parallel  to  the  Angular  Momentum  Axis 

In  this  mounting  scheme,  the  orientation  of  the  boresight  vector  is 
equal  to  the  orientation  of  the  angular  momentum  vector  as  illustrated  in 
Figure  5.7.  As  a  result,  estimation  of  the  boresight  vector  using  the  single¬ 
axis  SNR  approach  directly  estimates  the  angular  momentum  axis. 

This  is  the  preferred  mounting  scheme,  because  every  boresight  cal¬ 
culation  can  be  used  as  the  current  attitude  estimate.  Averaging  is  desir¬ 
able,  but  not  required.  No  assumptions  about  the  stability  are  needed. 
Estimates  can  be  continuously  made  throughout  torquing  maneuvers  that 
alter  the  position  of  the  angular  momentum  axis. 
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One  likely  scenario  is  that  the  angular  momentum  axis  is  relatively 
constant  over  a  short  time  period.  This  situation  occurs  between  torquing 
maneuvers.  Averaging  or  filtering  of  the  instantaneous  estimates  can  be 
employed  to  improve  the  accuracy.  Furthermore,  since  the  angular 
momentum  axis  is  constant  in  the  local  system,  its  orientation  can  be 
thought  of  in  terms  of  the  antenna  alignment  problem  as  describe  in  sec¬ 
tion  4.9.  In  this  case,  the  body  system  can  be  temporarily  defined  to  coin¬ 
cide  with  the  local  system  during  the  time  of  constant  orientation.  A  batch 
of  SNR  and  LOS  measurements  from  this  period,  can  be  used  in  Equation 
(4.21)  to  solve  for  the  elements  of  the  boresight  vector  in  the  body  frame, 
nB ,  where  we  set 

-  L  B 

H  =  n  . 

The  biggest  disadvantage  to  this  mounting  scheme  is  that  the  GPS 
antenna  has  to  be  placed  in  an  area  of  the  satellite  that  may  be  the  most 
desirable  for  other  instruments.  The  face  of  the  satellite  remains  in  an  ori- 
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entation  that  is  relatively  constant  compared  to  the  sides.  Scientific 
instruments  or  communications  antennas,  which  are  generally  the  pri¬ 
mary  mission  payload,  may  require  a  fixed  pointing  direction  and  thus, 
need  to  be  mounted  here.  A  lack  of  room  for  a  GPS  antenna  or  too  much 
signal  blockage  from  the  other  devices  might  be  the  result. 

5.5.2  Boresight  Mounted  Perpendicular  to  Angular  Momentum 
Axis 

In  this  mounting  scheme,  the  orientation  of  the  boresight  vector  sig¬ 
nificantly  changes  with  each  measurement  as  it  rotates  about  the  angular 
momentum  axis.  If  the  satellite  is  assumed  to  be  in  a  stable  spinning  ori¬ 
entation,  the  boresight  vector  always  remains  perpendicular  to  the  angular 
momentum  vector  as  shown  in  Figure  5.8.  Averaging  of  the  instantaneous 
boresight  estimates  is  not  applicable  because  there  is  no  component  of  the 
vectors,  excluding  the  errors  in  the  estimates,  that  lie  along  the  angular 
momentum  axis.  For  example,  suppose  local  system  is  defined  such  that 

(ltL)T  =  [o  0  i].  (5.23) 

The  boresight  vector  in  the  local  frame,  nL ,  is  going  to  be  of  the  form 

^  l  T  r  ~\ 

(nj )  =  coscotj  sincotj  0  , 

where  to  is  the  spin  rate.  For  perfect  and  equally  spaced  estimates,  averag- 
-  L 

ing  of  the  ni  values  results  in  the  null  vector.  For  non-perfect  estimates, 
the  average  still  tends  to  zero  if  the  estimation  errors  are  symmetric  about 
the  true  boresight.  Instead,  if  we  take  the  unit  vector  of  the  cross  product 
of  any  two  estimates  we  compute 
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L  T  l  T 

(AT)  x(Aj) 

0  0  coscotjsincotj- coscotjSincotj 

L  T  l  T 

(A“)  x(Aj) 

0  0  coscotjsincotj  -  coscotjSintotj 

(5.24) 


or 

(fiL)T=  [001].  (5.25) 

For  noisy  estimates,  we  must  take  the  average  over  a  batch  of  cross  prod¬ 
uct  results. 


This  mounting  scheme  has  several  disadvantages.  The  first  one  is 
that  batch  processing  of  the  instantaneous  estimates  is  a  necessity.  The 
angular  momentum  axis,  therefore,  must  remain  approximately  constant 
over  the  time  span  of  the  data  collection  batch.  This  condition  might  not  be 
met  during  maneuvers.  Other  problems  might  occur  if  there  is  too  much 
GPS  satellite  blockage  as  the  antenna  rotates  and  faces  the  direction  of  the 
Earth.  This  might  be  mitigated  by  weighting  the  boresight  estimates 
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according  to  their  satellite  visibility  conditions.  Estimates  made  with  the 
most  satellites  would  be  weighted  the  highest.  One  possible  approach  is 


y  F<Si’  SUk> 

i  =  1  n 


(5.26) 


where 


L  T  r  T| 

(Af)  XCA}*^) 


(5.27) 


L  T  l  T 

(hj )  and  (ni  +  k)  are  the  antenna  boresight  estimates  at  epoch  i  and  at 
an  epoch  k  measurements  later,  and  F(si,si+k)  is  some  empirically  deter¬ 
mined  function  of  the  number  of  satellites  in  view  at  the  two  epochs. 

In  late  1997,  the  SNOE  satellite,  described  earlier,  is  expected  to  be 
launched  aboard  a  Pegasus  XL  launch  vehicle.  It  has  a  GPS  antenna 
mounted  on  the  side  of  the  satellite  which  may  provide  data  to  analyze  this 
estimation  scenario. 


5.5.3  Boresight  Mounted  Acutely  to  the  Angular  Momentum  Axis 

In  this  mounting  scheme,  the  angle  between  the  boresight  vector 
and  the  angular  momentum  vector  is  acute  as  shown  in  Figure  5.9.  Since 
the  SNR  based  method  gives  only  single  axis  information,  knowing  the  ori¬ 
entation  of  the  antenna  boresight  in  the  body  frame  is  not  sufficient.  Aver¬ 
aging  or  some  type  of  batch  filtering  must  be  used.  As  the  antenna 
boresight  rotates  about  the  angular  momentum  axis,  it  maintains  a  con¬ 
stant  angle.  The  components  of  each  instantaneous  estimate  that  lie  in  the 
plane  perpendicular  to  the  angular  momentum  axis  average  out  to  zero  as 
the  boresight  cones  about  it.  The  components  that  lie  parallel  to  the  angu¬ 
lar  momentum  axis  average  to  a  constant  vector,  which  is  the  final  axis 
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estimate.  While  in  this  situation  the  boresight  orientation  in  the  local  sys¬ 
tem  is  not  constant,  the  average  of  equally  spaced  orientations  over  one 
revolution  is.  As  a  result,  the  boresight  alignment  technique  may  be  appli¬ 
cable. 


5.6  Three-Axis  SNR-Based  Attitude  Determination 

The  major  additional  requirement  for  the  three-axis  case  is  the  need 
for  either  another  antenna  or  measurable  azimuthal  variations  in  the 
antenna  gain  patterns.  Because  azimuthal  variations  are  very  difficult  to 
measure,  especially  at  the  level  of  accuracy  of  the  SNR  measurements,  this 
case  is  not  examined.  The  addition  of  a  second  non-parallel  antenna  adds 
the  estimation  of  another  axis  as  shown  in  Figure  5.10. 

L  B 

Computation  of  the  local-to-body  attitude  matrix,  C  ,  is  per¬ 
formed  in  two  steps.  The  first  involves  the  estimation  of  a  rotation  matrix, 
SCB ,  between  a  reference  frame,  denoted  by  the  superscript  S,  defined  by 
the  two  boresight  axes,  bi  and  b2,  and  the  body  frame,  denoted  with  the 
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Figure  5. 10.  Three  axis  estimation. 

superscript  B.  This  matrix  is  constant,  assuming  the  spacecraft  is  rigid 
and  the  body  position  of  the  antennas  do  not  change.  Using  the  alignment 
calibration  procedure  from  section  4.9,  the  two  boresight  vectors  are  esti¬ 
mated  in  the  body  frame.  An  orthonormal  set  of  axis  are  formed  to  arrive 
S  B 

at  C  as  described  in  section  2.5.  This  matrix  only  has  to  be  computed 
once. 

It  is  important  to  note  that  the  computation  for  depends  on  one 
of  three  conditions  being  satisfied.  Two  of  the  conditions  are  direct  require¬ 
ments  of  the  alignment  calibration  procedure.  As  in  the  single-axis  case, 
they  are  either  that  the  attitude  is  constant  or  that  a  truth  attitude  refer¬ 
ence  is  available  during  the  calibration  period.  Both  of  these  conditions 
are  readily  met  if  a  ground  calibration  is  performed.  If  the  mapping  func¬ 
tion  generated  from  the  antenna  gain  pattern,  then  the  third  condition  is 
that  the  precise  antenna  alignment  is  known  from  accurate  mechanical 
drawings. 


95 


The  second  step  involves  computing  the  transformation  between  the 

L  S 

local  system  and  the  boresight  system,  C.  The  single-axis  SNR-based 

estimation  computes  the  boresight  vectors  in  the  local  frame.  As  with  the 

first  transformation,  knowing  two  non-colinear  boresight  vectors  in  the 

L  S 

local  frame  allows  us  to  compute  an  orthonormal  set  of  axes  to  form  C. 

LB 

Combining  the  two  rotation  matrices  yields  C.  If  the  vehicle  attitude  is 
constant  over  the  data  collection  interval,  the  two  steps  can  be  combined  to 
directly  estimate  C  using  one  of  the  techniques  to  solve  for  the  minimi¬ 
zation  of  the  Wahba  cost  function  given  in  Equation  (2.22). 

5.7  Conclusions 

The  SNR-based  single-axis  estimation  technique  is  applicable  to 
both  the  estimation  of  the  angular  momentum  axis  of  a  spinning  satellite 
and  three-axis  estimation.  For  the  spinning  satellite  case,  if  the  antenna 
boresight  is  parallel  to  the  angular  momentum  axis,  the  SNR  method 
directly  estimates  it.  Assuming  the  axis  is  constant  allows  either  the  aver¬ 
aging  of  the  instantaneous  estimates  or  the  application  of  the  alignment 
calibration.  For  a  side  mounted  antenna,  the  cross  products  of  pairs  of  the 
instantaneous  estimates  have  to  be  averaged.  Variable  visibility  condi¬ 
tions  could  impose  the  need  for  weighting  of  the  estimates.  Three-axis 
estimation  requires  an  additional  antenna.  The  constant  transformation 
matrix  between  the  system  defined  by  the  antenna  boresights  and  body 
has  to  be  determined  from  either  a  calibration  set  or  from  mechanical 
drawings  of  the  antenna  orientations.  Each  pair  of  estimates  of  the  two 
antenna  boresights  are  used  to  estimate  the  transformation  matrix 
between  the  local  system  and  the  antenna  system. 


CHAPTER  6 


SNR  TEST  DATA 


6.1  Introduction 

The  data  used  in  this  research  for  SNR-based  single-axis  estimates 
are  derived  from  a  variety  of  receiver  and  antenna  combinations.  Each 
data  set  has  a  corresponding  source  of  truth  attitude  information  used  for 
performance  assessment  of  the  algorithms. 

Data  are  taken  from  five  different  sources.  Three  of  the  sources  are 
space  based  GPS  receivers  and  two  are  ground  receivers.  Three  of  the  data 
sources,  two  space  and  one  ground,  use  the  six  channel  Trimble  TANS  Vec¬ 
tor.  One  space  source  uses  the  eight  channel  Allen  Osborne  Associates  Tur- 
boStar  receiver.  The  other  ground  source  uses  a  12  channel  NovAtel 
receiver.  A  detailed  description  of  the  receivers  and  data  sets  follows. 

6.2  TANS  Vector  Receiver 

The  TANS  Vector  has  six  parallel  tracking  channels  and  a  multi¬ 
plexed  front  end  allowing  it  to  track  satellites  on  four  antennas  [Trimble, 
1993] .  Normally  with  four  antennas  in  operation,  it  measures  phase  differ¬ 
ences  between  one  antenna  designated  as  the  master  and  the  three  other 
antennas.  This  information,  taken  from  the  “E4”  data  packets,  along  with 
the  GPS  LOS  vectors,  taken  from  the  “EC”  data  packets,  are  used  to  esti- 
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mate  attitude.  In  this  paper,  measurements  from  only  one  antenna  are 
required.  A  picture  of  the  TANS  Vector  is  shown  in  Figure  6.1. 


antenna  inputs:  master,  #1,  #2,  #3 


Figure  6.1.  The  Trimble  TANS  Vector. 


6.3  NovAtel  12  Channel  Receiver 

The  NovAtel  receiver  is  based  on  the  NovAtel  GPSCard  OEM  Series 
receiver  model  number  3151R  [NovAtel,  1993].  The  GPSCard  Series  are 
multi-channel  parallel  tracking,  C/A  code  receivers.  Model  3151R  has  the 
capability  of  simultaneously  tracking  up  to  12  GPS  satellites.  The  antenna 
used  for  the  data  collection  in  this  research  is  a  microstrip  patch  antenna 
with  a  built-in  low  noise  amplifier  (LNA)  and  bandpass  filter.  Data  are  col¬ 
lected  using  a  PC  through  a  serial  port  connection.  For  this  research,  the 
“rnga”  data  packets  containing  the  SNR  measurements  in  dB  and  the 
“sata”  data  packets  containing  the  azimuth  and  elevation  to  the  GPS  satel- 
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lites  in  degrees  are  recorded.  A  picture  of  the  NovAtel  receiver  and 
antenna  is  shown  in  Figure  6.2. 


Figure  6.2.  The  12  Channel  NovAtel. 


6.4  GANE 

The  GANE  data  are  taken  from  the  GPS  Attitude  and  Navigation 
Experiment  (GANE),  flown  on  the  space  shuttle  Endeavour  during  May 
1996  [Carpenter  and  Hain,  1997].  The  measurements  are  from  a  TANS 
Vector  receiver  and  four  patch  antennas  mounted  on  a  structure  in  the 
shuttle  payload  bay.  The  antennas  are  mounted  in  a  rectangle  with  their 
boresights  pointing  out  of  the  payload  bay.  Each  antenna  is  surrounded  by 
a  choke  ring.  An  Inertial  Reference  Unit  (IRU)  provides  an  external  atti¬ 
tude  reference.  A  photograph  of  GANE  prior  to  its  placement  in  the  shut¬ 
tle  bay  is  shown  in  Figure  6.3. 
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Figure  6.3.  The  GANE  experiment.  [GANE  Home  Page] 

Data  are  from  two  different  orientations.  In  the  first  (GANE  #1), 
the  GANE  antennas  are  oriented  away  from  the  Earth  with  the  boresights 
in  the  zenith  direction  and  with  the  shuttle  approximately  tangent  to  its 
orbit.  In  the  second  (GANE  #2),  the  shuttle’s  nose  is  pitched  about  the 
orbit  normal.  The  GPS  antennas  are  oriented  approximately  41  deg  off 
zenith. 

Both  data  sets  analyzed  in  this  research  are  4  hours  long  with  SNR 
measurements  taken  approximately  every  5  seconds.  Each  set  contains 
about  a  2  hour  segment  when  a  reliable  truth  reference  is  available. 

6.5  CRISTA-SPAS 

The  CRISTA-SPAS  satellite,  built  by  Daimler-Benz  Aerospace,  was 
released  by  the  Space  Shuttle  Atlantis  in  November,  1994,  later  retrieved, 
and  returned  to  Earth.  The  primary  mission  was  to  explore  the  middle 
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atmosphere,  and  a  secondary  experiment  tested  GPS-based  attitude  deter¬ 
mination.  While  in  free  flight,  the  spacecraft  was  in  a  circular  orbit  at  an 
inclination  of  57  deg.  The  orbit  altitude  was  about  300  km.  The  satellite 
was  3-axis  controlled.  Primary  attitude  determination  was  provided  by  a 
star  tracker  and  gyro  inertial  measurement  unit  (IRU)  accurate  to  0.05  deg 
[Brock  et  al.,  1995]. 

In  the  secondary  experiment,  a  GPS  receiver  provided  by  Space  Sys¬ 
tems/Loral  was  used  to  independently  produce  attitude  data.  This  receiver 
was  a  TANS  Vector  modified  by  Stanford  University  for  attitude  determi¬ 
nation  in  space  [Lightsey  et  al.,  1994].  Four  GPS  patch  antennas  were 
mounted  on  the  zenith  face  of  the  spacecraft  [Ward,  L.M.,  1996]. 

The  total  data  set  from  CRISTA-SPAS  spans  a  period  of  over  30 
hours.  Measurements  are  at  intervals  of  3  to  10  seconds.  These  measure¬ 
ments  are  coordinated  with  the  IRU  attitude  and  with  precise  GPS  ephe- 
merides  computed  by  the  Jet  Propulsion  Laboratory  [Zumberge  and 
Bertiger,  1996] .  For  analysis,  the  IRU  data  is  used  as  a  truth  reference.  A 
picture  of  the  satellite  attached  to  the  space  shuttle  robot  arm  is  shown  in 
Figure  6.4. 

6.6  GPS/MET 

The  GPS-based  meteorology  (GPS/MET)  project  consists  of  an 
orbital  experiment  on  the  Orbital  Sciences  MicroLab-1  spacecraft.  The 
project  is  sponsored  by  the  University  Corporation  for  Atmospheric 
Research  (UCAR)  and  the  National  Science  Foundation  (NSF).  Launched 
on  April  3,  1995,  the  63  kg  satellite  operates  in  a  70  deg  inclined  orbit  and 
an  altitude  of  approximately  740  km.  Its  attitude  is  stabilized  by  a  gravity 


Figure  6.4.  The  CRISTA-SPAS  Satellite.  [JSC  Digital  Image 
Collection  Home] 

gradient  boom  and  is  controlled  by  torque  rods.  Attitude  estimates  are 
computed  from  the  filtering  of  data  from  Earth  sensors  and  a  magnetome¬ 
ter.  These  are  used  as  the  truth  reference.  Accuracy  is  expected  to  be 
between  5-10  deg  in  pitch  and  roll  and  10-15  deg  in  yaw  [Krebs,  1997]. 

On  board  is  a  dual  frequency  TurboStar  GPS  receiver  developed  by 
JPL  and  the  Allen  Osborne  Associates  (AOA).  Its  primary  objective  is  to 
take  high  frequency  measurements  of  signals  from  GPS  satellites  near  the 
Earth’s  horizon  to  sound  the  atmosphere.  Additionally,  1  Hz  measure¬ 
ments  are  taken  from  all  the  GPS  satellites  in  view.  A  diagram  of  the  sat¬ 
ellite  is  shown  in  Figure  6.5. 

For  this  work,  both  LI  and  L2  SNR  data  are  analyzed.  Three  data 
sets  collected  in  1997,  the  first  from  February  2  (GPS/MET  #1),  the  second 
from  February  20  (GPS/MET  #2),  and  the  third  from  February  26  (GPS/ 
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Figure  6.5.  Diagram  of  the  GPS/MET  satellite  [GPS/Met  Home  Page]. 

MET  #3),  are  examined.  Each  set  is  approximately  16  hours  in  duration. 
In  GPS/MET  #1  the  antenna  boresight  vector  is  approximately  in  the  sat¬ 
ellite  velocity  vector.  In  GPS/MET  #2  the  satellite  is  in  its  nominal  atti¬ 
tude  with  the  antenna  boresight  facing  opposite  the  velocity  vector.  In 
GPS/MET  #3,  the  antenna  boresight  is  in  the  orbit  normal  direction. 

6.7  Ground  TANS 

A  TANS  Vector  receiver  was  set  up  on  the  roof  of  the  University  of 
Colorado  Engineering  Center  as  shown  in  Figure  6.6.  It  was  attached  to 
four  antennas  mounted  on  a  rectangular  wooden  board.  The  orientation  of 
the  board  remained  constant  throughout  the  collection  period.  Data  was 
recorded  between  May  and  August,  1997. 

For  this  research,  data  sets  with  a  duration  of  at  least  20  hours  are 
analyzed.  Measurements  are  at  intervals  of  2.5  minutes.  The  SNR  data  is 
from  the  upper  right  antenna  which  is  designated  as  the  master.  The 
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phase-based  attitude  solution  from  all  four  antennas  is  used  as  the  truth 
reference.  Results  are  presented  from  a  set  referred  to  as  Ground  TANS. 


Figure  6.6.  Ground  experiment  setup. 


6.8  NovAtel  Ground 

A  12  channel  NovAtel  receiver  was  set  up  on  the  roof  of  the  Univer¬ 
sity  of  Colorado  Engineering  Center.  It  was  attached  to  the  center  of  a 
wooden  structure  as  shown  in  Figure  6.6.  As  with  the  Ground  TANS  sets, 
the  orientation  of  the  antenna  boresight  vector  was  constant.  Data  was 
collected  between  May  and  August,  1997. 

Data  sets  of  over  24  hours  in  duration  taken  at  intervals  of  2.5  min¬ 
utes  are  examined  in  this  research.  Data  from  three  different  antennas,  a 
NovAtel  patch,  a  Trimble  patch,  and  an  Ashtech  patch  with  a  choke  ring, 
are  used  for  the  various  data  collections.  Four  primary  data  sets  are  ana- 
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lyzed.  The  first,  called  NovAtel  #1,  comes  from  data  taken  with  the  NovA- 
tel  patch  antenna  in  May.  The  second,  called  NovAtel  #2,  comes  from  the 
Ashtech  antenna.  NovAtel  #3  comes  from  the  Trimble  antenna.  The  last 
set,  NovAtel  #4,  also  comes  from  the  NovAtel  antenna,  but  is  from  August. 
A  TANS  Vector  in  normal  attitude  mode  is  used  for  the  reference  attitude. 

6.9  a  vs.  SNR  for  the  Data  Sets 

Figures  6.7  through  6.13  show  the  a  versus  SNR  relationship  for  the 
various  data  sets.  For  the  purposes  of  a  comparison  between  all  of  the  data 
sets,  the  NovAtel  and  GPS/MET  SNR  output  are  converted  into  AMU’s 
using  Equation  (4.1).  Tb  get  the  Novatel  C/N0  to  the  same  range  as  the 
Vector  AMU’s,  a  bandwidth  of  333  Hz  is  assumed. 

To  quantify  the  overall  quality  of  the  relationship  between  off-bore- 
sight  angle  and  SNR,  the  value  of  Ap  in  degrees  is  computed.  It  is  defined 
as  the  mean  of  the  absolute  value  of  the  residuals  between  the  data  points 
and  a  5th  order  polynomial  fit  to  the  points. 
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Figure  6.13.NovAtel  a  vs.  SNR  (TANS  antenna). 


6.10  Normal  Distribution  Assumption 

Analysis  of  the  data  sets  shows  that  the  distribution  of  the  a’s  is 
symmetrically  distributed  in  each  SNR  bin.  The  distributions  vary  for 
each  data  set  and  for  each  SNR  bin.  Statistical  tests  indicate  that  the  dis¬ 
tributions  range  from  approximately  uniform  in  nature  to  Gaussian  or  in 
some  cases  a  contaminated  Gaussian  (a  Gaussian  with  a  specific  mean  and 
variance  mixed  with  a  Gaussian  having  a  different  mean  and  variance). 

6.10.1  Test  Statistic 

The  statistic  used  to  test  the  normal  distribution  assumption  is 
based  on  the  tail  weights  of  the  cumulative  distribution  function  (cdf) 
[Hoaglin,  et  al.,  1983].  The  analysis  shows  how  extreme  portions  of  the 
distribution  are  spread  out  in  relation  to  the  width  of  the  center.  The  test 
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statistic,  t,  is  defined  as  the  ratio  of  the  width  of  the  cdf  between  the 
median  and  the  99th  percentile  to  the  width  between  the  median  and  the 
75th  percentile.  This  ratio  is  then  normalized  by  the  ratio  for  a  true  nor¬ 
mal  distribution.  It  is  calculated  by 

I(F)  -F'1  (0.5)  /  (0.99)  -  O'1  (0.5)  (g  l) 

F"1  (0.75)  -  F'1  (0.5)  4T1  (0.75)  -4>_1  (0.5) 

where  F  is  the  cdf  of  the  test  data  and  4>  is  the  cdf  of  the  normal  distribu¬ 
tion.  Table  6.1  gives  the  values  of  x  for  several  common  distributions.  The 


Table  6.1.  x  for  some  standard  distributions.  [Hoaglin,  et  al.,  1983] 


Distribution 

x(F) 

Uniform 

0.57 

Triangular 

0.86 

Normal 

1.00 

CN(0.05;  3) 

1.20 

Double  -  Exponential 

1.63 

CN(0.05;  10) 

3.42 

Cauchy 

9.22 

distribution  CN(.05;  3)  represents  a  normal  distribution  with  j_l=0  and 
c2=l,  contaminated  by  a  5%  mixture  of  another  normal  distribution  with 
H=0  and  a2=3. 

Figures  6.14  through  6.18  show  a  plot  of  the  test  statistic  applied  to 
the  cdf  of  the  off-boresight  angles  for  each  SNR  bin  for  the  various  data 
sets.  Each  figure  consists  of  three  plots.  The  solid  line  corresponds  to  the 
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test  statistic  applied  to  each  SNR  bin  of  the  test  data.  The  dot-dash  line 
corresponds  to  the  test  statistic  applied  to  sets  of  normally  distributed  ran¬ 
dom  numbers  generated  in  MATLAB  with  mean,  variance,  and  set  size 
equal  to  those  of  the  test  data  bins.  The  dashed  line  corresponds  to  the  test 
statistic  applied  to  sets  of  uniformly  distributed  random  numbers  gener¬ 
ated  in  MATLAB  with  mean  and  set  size  equal  to  those  of  the  test  data 
bins  and  with  limits  equal  to  ±1.7  times  the  standard  deviation  of  the  test 
data  bins. 


Figure  6.14.  Distribution  test  for  GANE  #1. 
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While  the  data  clearly  show  some  contamination  from  a  true  normal 
distribution,  the  data  always  correspond  to  a  distribution  symmetric  about 
some  angle.  As  a  result,  methods  based  on  finding  the  maximum  likeli¬ 
hood  estimate  (MLE)  are  applicable. 

6.11  Conclusions 

The  variety  of  data  sets  allows  for  the  robustness  of  the  estimation 
algorithms  to  be  tested.  The  three  data  sets  taken  from  orbiting  satellites 
represent  the  most  realistic  data.  Fortunately,  these  sets  use  different  con¬ 
figurations  of  receivers  and  antenna  types  and  are  from  multiple  attitude 
orientations.  Data  taken  on  the  ground  allow  for  the  testing  of  a  wider 
range  of  receivers  and  antenna  combinations.  Furthermore,  the  collection 
of  data  using  the  TANS  Vector  on  the  ground,  permits  the  testing  of  a 
ground  calibration  set  in  conjunction  with  the  CRISTA-SPAS  orbit  data. 

Analysis  of  the  off-boresight  angle  versus  SNR  relationships  is  fur¬ 
ther  studied  in  the  next  chapter.  The  value  of  Ap  along  with  other  factors 
is  shown  to  be  a  good  indicator  of  the  estimation  performance. 

The  test  statistic,  x,  shows  that  the  a  versus  SNR  relationship  is 
approximately  symmetrically  distributed  in  the  SNR  bins  of  the  various 
data  sets.  In  the  next  chapter,  results  are  given  for  an  MLE  based  estima¬ 
tion  of  the  antenna  boresight  vector  involving  these  data  sets. 


CHAPTER  7 


SNR  BASED  SINGLE-AXIS  ESTIMATION:  MLE  RESULTS 

7.1  Introduction 

In  this  chapter  results  are  shown  for  the  estimation  of  a  single 
antenna  boresight  vector  using  the  MLE  algorithm  developed  in  Chapter  5 
applied  to  the  data  sets  described  in  the  previous  chapter.  The  accuracies 
range  from  3.0  to  10.5  deg  (RMS)  for  the  reference  results  and  4.4  to  11.1 
deg  (RMS)  for  results  using  various  mapping  function  sources.  These  val¬ 
ues  are  then  analyzed  according  to  a  number  of  performance  factors.  The 
most  important  are  shown  to  be  the  quality  of  the  mapping  function  and 
the  number  of  visible  measurements  per  epoch.  Based  on  the  performance 
demonstrated  from  all  of  the  data  sets,  this  method  is  predicted  to  have  a 
RMS  error  of  under  7  deg  for  a  high  quality  antenna  and  receiver  with  an 
average  visibility  of  5  or  more  satellites. 

The  best  results  are  from  the  GANE  satellite  data  which  had  a 
choke  ring  antenna.  The  choke  ring  contributes  to  a  good  mapping  func¬ 
tion.  The  NovAtel  results  benefit  by  the  12  channel  tracking  capability  of 
the  receiver  and  also  show  improvement  with  the  use  of  a  choke  ring 
antenna.  The  worst  results  are  from  the  GPS/MET  satellite.  Analysis  of 
these  results  is  hindered  by  a  poor  truth  reference. 
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7.2  Reference  Results 

Access  to  known  attitude  data  for  these  tests  allow  for  the  determi¬ 
nation  of  the  true  GPS  antenna  boresights.  First,  the  attitude  information 
is  used  to  rotate  the  LOS  vectors,  eL ,  recorded  in  the  local  frame  at  each 
epoch,  to  the  vehicle  body  frame,  eB .  These  vectors  are  collected  together 
for  a  batch  solution  of  the  boresight  alignment  in  the  body  frame,  nB ,  as 
detailed  in  Section  4.9.  Next,  the  attitude  data  are  used  to  rotate  nB  to  nL 
at  each  epoch.  Finally,  the  corresponding  a  values  are  calculated  from  the 
angles  between  nL  and  eL  from  which  a  mapping  function  is  generated. 
The  mapping  functions  for  each  satellite  and  antenna  combination  are 
shown  in  Figures  7.1  through  7.6. 

The  computed  alignments  for  each  data  set  are  shown  in  Table  7.1.  It 
is  expected  that  the  estimated  alignments  between  the  space-based  data 
sets  for  each  receiver  should  be  the  same.  The  ground  sets  are  not  expected 
to  be  equal  to  each  other  because  the  mounting  conditions  are  different  for 
each  antenna.  The  alignments  for  the  two  CRISTA-SPAS  sets  are  within 

2.3  deg  of  each  other.  The  GPS/MET  alignments  are  all  within  2.6  deg.  The 
two  GANE  alignments,  however,  are  over  10  deg  apart.  After  combining  the 
two  data  sets,  the  computed  alignment  agrees  to  within  0.7  deg  of  the  bore- 
sight  estimated  with  just  the  GANE  #1  set.  Considering  the  fact  that  the 
GANE  #2  set  covers  a  larger  range  of  angles  with  fewer  data  points  than 
the  GANE  #1  set,  the  results  suggest  that  there  may  not  be  sufficient  data 
for  an  accurate  alignment  calibration  from  just  the  GANE  #2  set. 


a  (deg)  a  (deg) 
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Figure  7.1  GANE  mapping  function. 


Figure  7.2  CRISTA-SPAS  mapping  function. 
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Figure  7.4  NovAtel  mapping  function  (NovAtel  antenna) 
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Figure  7.5  NovAtel  mapping  function  (Ashtech  antenna). 


Figure  7.6  Ground  TANS  mapping  function. 
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Table  7.1.  Boresight  Alignment. 


Data  Set 

boresight  coordinates  in  the  body  frame,  nB 

X 

y 

z 

GANE  #1 

0.0113 

-0.0213 

0.9997 

GANE  #2 

-0.0406 

0.1496 

0.9879 

GANE  #1  and  #2 

0.0157 

-0.0093 

0.9998 

CRISTA-SPAS  #1 

-0.9940 

0.0810 

-0.0731 

CRISTA-SPAS  #2 

-0.9939 

0.0505 

-0.0984 

GPS/MET  #1 

0.9964 

-0.0136 

-0.0842 

GPS/MET  #2 

0.9915 

-0.0277 

-0.1272 

GPS/MET  #3 

0.9949 

-0.0152 

-0.0996 

Ground  TANS 

0.0282 

0.0129 

0.9995 

NovAtel  #1  (NovAtel) 

-0.0062 

-0.0658 

0.9978 

NovAtel  #2  (Ashtech) 

0.0051 

-0.0095 

0.9999 

NovAtel  #3  (TANS) 

0.0636 

-0.0166 

0.9978 

NovAtel  #4  (NovAtel) 

-0.0362 

-0.0188 

0.9992 

The  MLE  solution,  as  described  in  Chapter  5,  is  used  to  estimate  a 
new  n  at  each  measurement  epoch.  Table  7.2  shows  the  RMS  error,  maxi¬ 
mum  error,  and  the  90th  percentile  error  for  each  data  set  described  in  the 
last  chapter.  For  these  results  the  same  data  used  to  create  the  reference 
calibration  are  used  to  estimate  n .  Figures  7.7  through  7.10  show  the  esti¬ 
mation  errors  as  a  function  of  time  for  four  of  the  data  sets,  one  each  from 
GANE,  CRISTA-SPAS,  GPS/MET,  and  the  NovAtel  ground  sets.  Each  solu¬ 
tion  is  based  upon  a  single  measurement  epoch. 
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Table  7.2.  Estimation  Results  with  Reference. 


Data  Set 

RMS  Error  (deg) 

90%  Error  (deg) 

Max.  Error  (deg) 

GANE  #1 

3.0 

4.7 

8.6 

GANE  #2 

5.9 

8.5 

23.9 

CRISTA-SPAS  #1 

6.4 

9.5 

23.1 

CRISTA-SPAS  #2 

7.1 

10.8 

21.2 

GPS/MET  #1 

10.5 

15.9 

43.9 

GPS/MET  #2 

10.4 

15.8 

52.6 

GPS/MET  #3 

10.0 

15.4 

52.6 

Ground  TANS 

4.8 

7.7 

15.2 

NovAtel  #1 

6.1 

9.4 

31.1 

NovAtel  #2 

4.3 

6.4 

30.2 

NovAtel  #3 

6.0 

8.9 

45.0 

NovAtel  #4 

6.1 

9.2 

20.8 

The  results  for  GANE  #1  are  clearly  the  best.  Not  only  is  the  RMS 
error  under  4  deg  but  the  maximum  error  is  under  9  deg.  NovAtel  #2,  which 
uses  an  Ashtech  choke  ring  antenna,  is  the  second  best  data  set  in  terms  of 
RMS  error,  but  it  has  a  large  maximum  error.  The  90th  percentile  errors  fol¬ 
low  the  same  trend  as  the  RMS  errors.  Recall  from  Chapter  5  that  the 
NovAtel  #1  and  #3  sets  use  a  NovAtel  patch  and  Trimble  patch  antenna, 
respectively.  A  comparison  between  their  results  shows  no  significant  dif¬ 
ference  between  their  performances.  The  Ground  TANS  set,  which  is  taken 
in  exactly  the  same  environment  as  the  NovAtel  sets  and  uses  the  Trimble 
patch  antenna,  produces  results  that  are  almost  as  good  as  the  NovAtel  set 
taken  with  the  choke  ring  antenna.  All  of  the  data  sets,  with  the  exception  of 
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GANE  #1,  exhibit  spikes  in  the  error  plots.  The  maximum  errors  occur  at 
these  spikes.  The  largest  spikes,  as  evident  in  both  the  GPS/MET  and 
NovAtel  plots,  are  significantly  larger  than  the  general  trend  of  the  errors. 
Since  the  NovAtel  data  is  taken  on  the  ground  and  consists  of  several  days 
worth  of  measurements,  a  repetition  of  the  test  conditions  occurs.  The 
spikes  are  shown  to  be  repeated  at  the  time  of  each  day  when  the  satellite 
geometry  is  identical.  Further  analysis  for  the  causes  of  these  bad  points  is 
given  in  Section  7.5  where  the  performance  factors  are  discussed. 

In  Table  7.3  the  estimate  of  n  is  equal  to  the  initial  estimate,  nQ,  of 
the  gradient  search  as  given  by  Equation  (5.10).  The  RMS  errors  for  most 
of  the  data  sets  are  within  0.3  deg  of  the  Table  7.2  results.  The  GANE  #2 
results,  however,  are  significantly  larger.  The  GANE  #2  set  has  the  lowest 
average  number  of  satellites  in  view  with  3.8  compared  to  over  5.4  for  the 
other  sets.  This  suggests  that  the  gradient  search  is  only  useful  when  the 
number  of  measurements  is  small.  For  a  large  number  of  measurements, 
the  normalized  result  of  the  weighted  least  squares  solution  is  sufficient. 
The  GPS/MET  results,  on  the  other  hand,  are  up  to  0 .9  deg  smaller.  This  dis¬ 
crepancy  may  be  caused  by  the  poor  truth  reference. 

7.3  Results  for  Specific  Mapping  Function  Sources 

In  Table  7.4  the  data  used  to  estimate  n  are  independent  of  the  data 
used  to  generate  the  calibration  function.  The  GANE  data  are  evaluated 
with  both  a  gain  pattern  mapping  function  and  by  using  an  alternate  set  of 
calibration  data.  The  CRISTA-SPAS  data  are  evaluated  with  a  gain  pat¬ 
tern  function,  a  ground  based  calibration  function,  and  an  alternate  space 
based  generated  function.  The  GPS/MET  data  is  tested  on  calibration 
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Table  7.3.  Estimation  Results  with  Reference  using  n0 . 


Data  Set 

RMS  Error  (deg) 

90%  Error  (deg) 

Max.  Error  (deg) 

GANE  #1 

3.2 

4.7 

21.4 

GANE  #2 

8.7 

8.7 

71.6 

CRISTA-SPAS  #1 

6.2 

9.5 

20.5 

CRISTA-SPAS  #2 

7.1 

10.8 

17.0 

GPS/MET  #1 

9.9 

15.1 

32.7 

GPS/MET  #2 

9.3 

14.1 

52.5 

GPS/MET  #3 

9.7 

14.9 

48.8 

Ground  TANS 

5.0 

7.5 

22.7 

NovAtel  #1 

6.1 

9.4 

25.1 

NovAtel  #2 

4.2 

6.4 

21.5 

NovAtel  #3 

5.7 

8.6 

33.8 

NovAtel  #4 

6.0 

9.1 

18.5 

functions  created  from  alternate  data.  The  NovAtel  data  is  used  to  test  the 
variability  of  the  mapping  functions  over  an  extended  time  period.  To 
accomplish  this,  the  mapping  function  generated  from  NovAtel  #1  data, 
collected  over  two  months  before  NovAtel  #4,  is  used  in  conjunction  with 
this  later  set.  Listed  to  the  right  of  the  error  values  are  the  differences 
from  the  results  in  Table  7.2. 

In  all  of  the  cases,  both  the  RMS  errors  and  the  90th  percentile 
errors  increase  over  the  reference  results.  However,  the  absolute 
increases,  are  not  very  large.  They  ranged  from  0.1  deg  to  2.0  deg  for  the 
RMS  errors  and  0.6  deg  to  2.9  deg  for  the  90th  percentile  errors.  The 
GANE  #1  data  set  is  the  most  sensitive  to  the  change  in  the  mapping  func- 
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Table  7.4.  Estimation  Results  with  Various  Mapping  Functions. 


Data  Set 

RMS  Error  (A) 
(deg) 

90%  Error  (A) 
(deg) 

Max.  Error  (A) 
(deg) 

GANE  #1  w/  gain 
pattern 

4.6  (+1.6) 

7.3  (+2.6) 

11.7  (+3.1) 

GANE  #1  w/ 

GANE  #2  data 

5.0  (+2.0) 

7.6  (+2.9) 

17.6  (+9.0) 

GANE  #2  w/  gain 
pattern 

7.4  (+1.5) 

10.8  (+2.3) 

25.4  (+1.5) 

CRISTA-SPAS  #1 
w/  gain  pattern 

6.6  (+0.2) 

9.7  (+0.2) 

27.4  (+4.3) 

CRISTA-SPAS  #1 
w/  ground  data 

7.2  (+0.8) 

10.8  (+1.3) 

24.3  (+1.2) 

.CRISTA-SPAS  #2 
w/CRISTA- 
-SPAS  #1  data 

7.7  (+0.6) 

11.6  (+0.8) 

22.3  (+1.1) 

GPS/MET  #3  w/ 
GPS/MET  #1  data 

10.7  (+0.7) 

10.7  (+0.7) 

63.2  (+10.6) 

GPS/MET  #lw/ 
GPS/MET  #3  data 

11.1  (+0.6) 

16.6  (+0.7) 

42.0  (-1.9) 

Ground  TANS  w/ 
CRISTA-SPAS  #1 
data 

5.9  (+1.1) 

9.5  (+1.8) 

17.2  (+2.0) 

NovAtel  #4  w/ 
NovAtel  #1  data 

7.0  (+0.9) 

10.9  (+1.4) 

32.0  (12.8) 

NovAtel  #2  w / 
gain  pattern 

4.4  (+0.1) 

6.5  (+0.1) 

22.2  (-8.0) 

tion,  while  the  GPS/MET  sets  are  the  least  sensitive.  This  is  possibly  due 
to  the  GANE  #1  data  having  the  best  data  quality,  and  as  a  result,  is  more 
sensitive  to  small  changes  in  the  calibration  function.  GPS/MET,  on  the 
other  hand,  has  the  worst  data  quality  and  as  a  result,  is  not  sensitive  to 
small  changes  in  the  calibration  function.  The  use  of  gain  pattern  based 
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mapping  functions  works  as  well  as  the  use  of  alternate  calibration  data. 
The  maximum  errors  fluctuated  unpredictably  from  a  drop  of  more  than 
15  deg  to  an  increase  of  over  12  deg. 

7.4  CRISTA-SPAS  Results  using  Individual  Satellite  Mapping 
Functions 

Analysis  of  the  plots  presented  in  Section  4.7  shows  that  there  is  no 
significant  difference  between  a  mapping  function  based  on  all  of  the  GPS 
satellites  and  functions  based  on  individual  ones.  Table  7.4  presents  results 
for  the  CRISTA-SPAS  satellite  demonstrating  this.  The  individual  mapping 
functions  are  generated  from  the  CRISTA-SPAS  #1  set.  When  compared  to 
the  single  mapping  function  reference  results,  the  reference  results  based  on 
individual  functions  are  significantly  better.  This  is  expected  because  the 
individual  satellites  functions  are  based  on  a  much  smaller  range  of  data.  As 
a  result,  each  SNR  and  a  pair  are  very  closely  correlated  with  the  average 
used  in  the  specific  function.  However,  this  is  not  the  case  for  the  values  in 
the  CRISTA-SPAS  #2  set.  The  results  when  using  the  single  mapping  func¬ 
tion  from  the  CRISTA-SPAS  #1  set  are  clearly  superior. 

7.5  Performance  Factors 

The  results  shown  in  the  previous  sections  illustrate  the  effect  of 
mapping  function  errors  on  the  achievable  accuracy.  Other  factors  which 
affect  the  performance  are  the  quality  of  SNR  data  as  measured  by  its  Ap 
value  and  the  choice  of  SNR  bin  size.  These  factors  are  governed  by  the 
data  processing  technique.  Additionally,  the  performance  is  affected  by  the 
total  number  of  satellites  tracked  and  the  computed  PDOP  value  which  are 
both  independent  of  the  processing  technique. 
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Table  7.5.  Estimation  Results  with  Individual  Mapping  Functions. 


Data  Set 

RMS  Error  (A) 
(deg) 

90%  Error  (A) 
(deg) 

Max.  Error  (A) 
(deg) 

CRISTA-SPAS  #1 

6.4 

9.5 

23.1 

w/  single  mapping 
function 

(reference  results) 

CRISTA-SPAS  #1 

w/  individual  satellite 
mapping  functions 

(reference  results) 

4.7 

7.3 

18.8 

CRISTA-SPAS  #2 

7.7 

11.6 

22.3 

w/  single  mapping 
function  from  CRISTA- 
SPAS  #1 

CRISTA-SPAS  #2 

11.3 

16.6 

73.3 

w/  individual  satellite 
mapping  functions  from 
CRISTA-SPAS  #1 

Table  7.6  lists  the  quality  statistics  for  each  data  set.  A  comparison 
of  the  data  quality  statistics  and  the  estimation  results  in  Tables  7.2  and  7.4 
yields  some  interesting  insights.  The  two  GANE  data  sets  have  similar 
measurement  error  distributions  as  evidenced  by  the  Ap  values,  but  signif¬ 
icantly  different  satellite  visibilities  (mean  of  5.4  to  3.8  satellites).  The 
resulting  RMS  attitude  errors  are  3.0  and  6.4  deg  for  GANE  #1  and 
GANE  #2,  respectively.  CRISTA-SPAS  #1  has  about  the  same  mean  num¬ 
ber  of  satellites  as  GANE  #1,  but  the  measurement  error  factor  Ap  is  7. 1  as 
compared  to  3.8  for  GANE  #1.  The  resulting  RMS  attitude  error  for 
CRISTA-SPAS  #1  is  6.4  deg.  The  NovAtel  sets  have  relatively  poor  map¬ 
ping  functions.  However,  the  ability  of  the  receiver  to  track  up  to  12  satel¬ 
lites  greatly  enhances  the  satellite  visibility  and  consequently  the  attitude 
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Table  7.6. 

Quality  of  Data 

Data  Set 

Ap  (deg) 

Mean  #  of  Sats 

Mean  PDOP 

GANE  #1 

3.8 

5.4 

2.3 

GANE  #2 

4.3 

3.8 

32.0 

CRISTA-SPAS  #1 

7.1 

5.8 

2.2 

CRISTA-SPAS  #2 

7.3 

5.6 

2.4 

GPS/MET  #1 

9.4 

7.4 

1.5 

GPS/MET  #2 

10.4 

7.3 

1.8 

GPS/MET  #3 

9.3 

6.9 

2.1 

Ground  TANS 

5.8 

5.8 

2.2 

NovAtel  #1 

10.0 

8.2 

1.2 

NovAtel  #2 

7.7 

8.2 

1.2 

NovAtel  #3 

10.0 

8.2 

1.3 

NovAtel  #4 

9.9 

7.9 

1.3 

estimation  performance.  The  GPS/MET  data  have  quality  that  is  similar 
to  the  NovAtel  sets.  The  satellite  visibility,  on  the  other  hand,  averages  to 
about  1  fewer  satellites  per  measurement  epoch.  As  a  result,  the  conse¬ 
quences  of  poor  data  quality  are  not  mitigated  as  well. 
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Tables  7.7  through  7.10  show  a  break  down  of  the  estimation  errors 
versus  the  number  of  satellites  in  view  for  the  combined  GANE  data  sets, 
CRISTA-SPAS  #1,  GPS/MET  #1,  and  NovAtel  #1,  respectively.  To  generate 
the  results,  solutions  are  formed  using  only  the  specified  number  of  satel¬ 
lites.  The  first  m  satellites  are  selected  rather  than  choosing  them  based 
on  a  geometrical  constraint.  This  is  repeated  over  the  possible  range  of 
satellite  visibilities. 


Table  7.7.  Visibility  Dependence  for  GANE  #1  and  #2 


Data  Statistic 

Number  of  Satellites 

3 

4 

5 

6 

RMS  Error  (deg) 

7.0 

3.3 

3.2 

2.9 

90%  Error  (deg) 

6.8 

5.3 

4.8 

4.4 

Max.  Error  (deg) 

68.4 

10.4 

8.6 

8.0 

Table  7.8.  Visibility  Dependence  for  CRISTA-SPAS  #1 

Data  Statistic 

Number  of  Satellites 

3 

4 

5 

6 

RMS  Error  (deg) 

13.8 

10.4 

6.9 

5.6 

90%  Error  (deg) 

20.8 

13.9 

10.5 

8.3 

Max.  Error  (deg) 

70.4 

84.4 

23.8 

20.6 

In  every  data  set  and  for  each  one  of  the  computed  statistics,  the 
errors  decrease  as  more  satellites  are  used.  The  largest  decrease  over  the 
span  of  visibility  conditions  occur  with  the  data  sets  having  the  worst  data 
quality.  The  largest  decreases  occur  in  the  GANE  and  CRISTA-SPAS  sets 
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Table  7.9.  Visibility  Dependence  for  GPS/MET  #1 


Data  Statistic 

Number  of  Satellites 

5 

6 

7 

8 

RMS  Error  (deg) 

11.9 

11.2 

10.3 

9.6 

90%  Error  (deg) 

18.4 

17.1 

16.0 

14.8 

Max.  Error  (deg) 

42.0 

43.9 

37.1 

32.5 

Table  7.10.  Visibility  Dependence  for  NovAtel  #1 


Number  of  Satellites 

UCLxjCL  UtatlotlV/  - 

5 

6 

7 

8 

9 

10 

11 

RMS  Error  (deg) 

9.0 

8.4 

7.5 

6.5 

6.2 

6.0 

4.5 

90%  Error  (deg) 

13.3 

12.7 

11.4 

10.0 

9.6 

8.7 

3.8 

Max.  Error  (deg) 

57.9 

38.4 

33.5 

19.0 

16.7 

13.0 

5.9 

when  the  visibility  is  increased  from  3  to  4  satellites.  The  visibility  depen¬ 
dences  clearly  show  the  advantage  of  a  multichannel  receiver  with  the 
capability  of  tracking  a  large  number  of  satellites,  especially  when  the 
quality  of  the  mapping  function  is  poor. 

Figures  7.11  through  7.14  are  scatter  plots  showing  the  relationship 
between  the  estimation  error  and  PDOP  value.  In  Figure  7.11,  the  PDOP’s 
that  are  greater  than  16  are  grouped  together  for  better  resolution  in  the 
graph.  None  of  the  plots  demonstrate  a  strong  correlation  between  these 
two  values.  The  errors  from  the  GANE  data  are  all  under  20  deg  except 
when  the  PDOP  is  larger  than  16.  The  CRISTA-SPAS  data  shows  a  small 
increase  in  the  maximum  errors  when  the  PDOP  is  over  3.  The  largest 
NovAtel  errors  occur  for  PDOP’s  greater  than  1.7.  GPS/MET  has  the  least 
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PDOP  dependence  of  the  analyzed  data  sets.  In  all  cases,  however,  a  large 
PDOP  does  not  guarantee  a  large  estimation  error.  Each  data  sets  demon¬ 
strates  instances  of  low  errors  with  a  much  higher  than  average  PDOP. 
While  it  might  be  expected  that  large  PDOFs,  indicating  poor  geometry, 
would  be  a  major  contributor  to  poor  estimation,  the  data  does  not  demon¬ 
strate  this. 

Figures  7.15  and  7.16  show  the  dependence  on  SNR  bin  resolution 
for  the  NovAtel  and  TANS  Vector  receivers,  respectively.  For  both  receiv¬ 
ers  the  RMS  errors  show  a  slight  but  steady  increase  as  the  bin  size  gets 
larger.  The  maximum  errors  generally  increase  with  bin  size  but  there  is 
some  fluctuation.  In  both  cases,  an  extremely  small  bin  size  is  not  dramat¬ 
ically  beneficial,  but  does  show  some  improvement.  It  was  thought  that 
SNR  variations  produced  by  multipath  would  indicate  a  larger  bin  size; 
however,  experimental  results  do  not  bear  this  out.  Since  the  mapping 
functions  are  generated  prior  to  estimation,  there  is  no  reason  on  the  basis 
of  computational  speed  to  use  larger  bin  sizes. 

7.6  Estimation  Performance  vs.  Computed  Error  Bounds 

In  Section  5.4.3  Equations  (5.21)  and  (5.22)  are  derived,  giving  a 
means  by  which  the  expected  errors  for  the  two  boresight  angles  are  deter¬ 
mined  at  each  epoch.  These  values  correspond  to  a  prediction  of  the  actual 
error  to  a  la  level  of  certainty.  In  other  words,  it  is  expected  that  about 
67%  of  the  angle  errors,  80  and  5<J>,  computed  for  the  reference  results  will 
be  less  than  or  equal  to  these  predicted  values  at  a  given  epoch. 

Table  7.11  shows  a  comparison  of  the  actual  angle  errors  to  the  pre¬ 
dicted  errors.  The  predicted  error  bounds  generally  agree  to  within  la  of 


Max  Error  (deg)  RMS  Error  (deg)  era  Max  Error  (deg)  RMS  Error  (deg) 
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o-o  GANE  #1  x-x  Crista-Spas  #1 


ire  7.15  SNR  bin  resolution  dependence  for  GANE  and 
CRISTA-SPAS  #1. 


o-o  NovAtel  Patch  x-x  NovAtel  Choke 


Figure  7.16  SNR  bin  resolution  dependence  for  NovAtel  #1  and  #2 
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Table  7.11.  Comparison  of  Actual  Errors  to  Predicted  Errors 


Data  Set 

RMS  Error 
(deg) 

%  of  50’s  within 
predicted  bounds 

%  of  8<t>’s  within 
predicted  bounds 

GANE  #1 

3.0 

67 

76 

GANE  #2 

6.4 

66 

63 

CRISTA-SPAS  #1 

6.4 

69 

72 

CRISTA-SPAS  #2 

7.1 

-  60 

70 

GPS/MET  #1 

10.5 

59 

49 

GPS/MET  #2 

10.4 

57 

58 

GPS/MET  #3 

10.0 

59 

56 

Ground  TANS 

4.8 

64 

63 

NovAtel  #1 

6.1 

73 

63 

NovAtel  #2 

4.3 

68 

73 

NovAtel  #3 

5.9 

68 

73 

NovAtel  #4 

6.1 

72 

71 

the  actual  errors.  The  data  sets  with  the  smallest  errors,  GANE  #1, 
CRISTA-SPAS  #1,  and  NovAtel  #2,  show  the  best  agreement.  The  biggest 
exception  is  with  the  GPS/MET  data  which  shows  a  much  smaller  percent¬ 
age  of  estimates,  as  low  as  49%,  within  the  bounds.  One  possible  reason 
for  this  discrepancy  is  that  the  GPS/MET  data  is  analyzed  using  a  poor 
truth  reference,  implying  that  the  mapping  function  errors  are  likely  to  be 
correlated  or  non-Gaussian. 

7.7  Conclusions 

This  chapter  has  presented  the  results  for  an  MLE  based  approach 
to  estimate  a  single  axis  of  attitude  using  the  SNR  measurements  from  a 
single  GPS  antenna.  This  method  is  based  upon  an  empirical  mapping  of 
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the  angle  between  the  GPS  antenna  off-boresight  angle  and  the  SNR  mea¬ 
surement  of  the  GPS  signal.  The  RMS  of  the  estimation  errors  range  from 
3.0  to  10.5  deg,  the  90th  percentile  errors  range  from  4.7  to  15.9  deg,  and 
the  maximum  errors  range  between  8.6  and  52.6  deg.  The  best  results 
occur  when  there  is  a  good  mapping  function  quantified  by  the  mean  resid¬ 
ual  value  and  when  there  are  6  or  more  satellites  in  view.  The  best  results 
are  from  GANE  #1  that  uses  an  antenna  and  choke  ring  combination  and 
is  in  an  orientation  that  allows  for  high  visibility.  The  second  best  results 
are  from  NovAtel  #2  that  also  uses  an  antenna  and  choke  ring  combina¬ 
tion.  While  the  mean  residual  value  of  the  NovAtel  mapping  functions  are 
not  as  good  as  the  ones  for  the  GANE  data  sets,  the  receiver’s  ability  to 
track  up  to  12  satellites  at  once  allows  for  a  significantly  greater  mean 
number  of  satellites  in  view.  The  worst  results  are  from  the  GPS/MET  sat¬ 
ellite.  The  data  sets  have  a  relatively  good  visibility,  but  their  mapping 
functions  are  among  the  worst  of  all  the  sets.  Analysis  of  these  sets  is  hin¬ 
dered  by  a  lack  of  a  good  truth  reference.  Without  a  good  reference,  the 
calibration  functions  are  degraded  and  the  computed  errors  are  not  accu¬ 
rate. 

As  expected,  the  estimation  performance  is  degraded  when  the 
source  of  calibration  data  is  not  the  actual  data  being  used.  However,  the 
largest  increase  in  the  RMS  error  is  only  2  deg  and  the  largest  increase  in 
the  90th  percentile  error  is  2.9.  The  gain  pattern  based  mapping  functions 
work  as  well  as  calibration  data  based  ones. 

The  analysis  of  the  performance  factors  clearly  show  that  the  qual¬ 
ity  of  the  mapping  function  as  given  by  its  Ap  value  and  the  number  of  sat- 
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ellites  being  used  are  the  two  best  indicators  of  performance.  There  is  only 
a  slight  indication  that  PDOP  has  a  noticeable  effect. 

The  results  demonstrated  in  this  research  show  that  this  method 
has  the  potential  for  coarse  single  attitude  estimation  to  a  level  of  5  to 
10  deg  with  one  antenna.  The  use  of  various  calibration  sources  and  the 
satisfactory  estimation  results  indicate  that  the  technique  can  be  practi¬ 
cally  implemented.  In  particular  it  would  be  useful  in  the  role  of  emer¬ 
gency  backup  or  integer  ambiguity  resolution. 


CHAPTER  8 


SUMMARY 


Two  sets  of  techniques  have  been  presented  in  this  research  that  use 
GPS  information  for  spacecraft  attitude  estimation.  In  particular,  these 
methods  are  designed  to  meet  the  criteria  of  low  cost  satellite  missions. 
The  first  set  of  methods,  based  on  GPS  carrier  phase  data,  can  only  be 
applied  to  satellites  using  spin  stabilization.  The  second  set  of  methods, 
based  on  GPS  SNR  data,  can  be  used  for  any  type  of  satellite  stabilization. 

The  carrier  phase  approach  is  based  on  the  phase  differences 
between  two  antenna  mounted  on  one  face  of  a  spinning  satellite.  Because 
the  antennas  experience  large  changes  in  position  in  a  short  time  span,  the 
measurements  can  be  differenced  in  time  for  the  computation  of  a  displace¬ 
ment  vector  observable.  This  calculation  is  not  significantly  affected  by  the 
carrier  phase  integer  ambiguity.  As  a  result,  it  is  not  computationally 
intensive  and  is  not  hindered  by  cycle  slips. 

Both  the  satellite  frequency  and  attitude  information  are  extracted 
from  the  displacement  vectors.  The  frequencies  are  determined  by  the  AR 
and  FFT  frequency  domain  methods  and  a  time  domain  technique  that  is 
based  on  the  average  of  the  displacement  vectors.  The  AR  method  is  the 
most  accurate  of  the  three.  For  a  batch  of  measurements  spanning  a 
minute,  it  produces  estimates  to  within  0.03  rpm  for  a  5  rpm  satellite. 
After  five  minutes  the  error  is  under  0.01  rpm.  The  orientation  of  the  sat- 
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ellite  angular  momentum  vector  is  estimated  by  a  batch  filter  and  an  aver¬ 
aging  technique  using  the  cross  products  of  the  displacement  vectors.  Both 
the  cross  product  and  batch  filter  techniques  produce  pointing  errors  of 
under  0.25  deg.  In  addition,  the  batch  filter  estimates  the  nutation  angle 
to  within  0.2  deg  and  the  phase  angles  to  within  2  deg. 

The  SNR-based  approaches  estimate  the  antenna  boresight  vector 
by  using  the  relationship  between  the  SNR  of  a  GPS  measurement  and  the 
corresponding  angle  a.  The  a’s  are  between  the  boresight  and  the  LOS 
vectors.  Because  of  the  variations  in  the  antenna  gain  patterns,  the  SNR 
is  highest  at  small  a’s  and  decreases  as  a  increases.  A  techniques  has  been 
presented  using  this  relationship.  It  computes  the  boresight  vector  using  a 
maximum  likelihood  estimation  (MLE). 

The  MLE  method  requires  the  creation  of  a  mapping  function 
between  SNR  measurements  and  the  corresponding  a’s.  This  is  formed  by 
either  using  a  calibration  data  set  from  SNR’s  taken  at  known  a’s  or  by 
using  a  gain  pattern  measurement  scaled  by  the  receiver  SNR  limits.  The 
calibration  approach  has  the  advantage  of  using  data  that  is  more  similar 
to  the  actual  estimation  data,  while  the  gain  pattern  approach  has  the 
advantage  of  requiring  less  initialization  data.  After  the  mapping  func¬ 
tions  are  created,  axis  estimates  at  each  epoch  are  formed  using  the  LOS 
vectors,  the  a  values  mapped  from  each  of  the  SNR  measurements,  and  the 
weights  assigned  to  the  a  values. 

The  SNR  vs.  a  relationship  is  affected  by  a  number  of  influences. 
Two  of  these,  the  space  loss  and  gain  pattern  of  the  GPS  transmitting 
antenna,  are  adjusted  for  by  knowing  the  satellite  positions.  Power  fluctu¬ 
ations  among  GPS  satellites  do  not  produce  a  significant  effect.  Multipath 
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causes  the  largest  fluctuations  in  the  relationship.  This  effect  is  treated  as 
measurement  noise. 

The  MLE  method  has  been  applied  to  a  variety  of  data,  including 
both  ground  and  space  sources.  The  best  result  of  3.0  deg  (RMS)  is  from 
the  GANE  satellite.  The  worst  result  of  10.5  deg  (RMS)  is  from  the  GPS/ 
MET  satellite.  The  results  show  that  two  factors  contribute  most  signifi¬ 
cantly  to  the  performance  of  the  technique.  The  first  is  the  quality  of  the 
SNR  data.  A  good  quality  consists  of  a  small  potential  range  of  oc’s  for  each 
SNR  value.  Choke  ring  antennas  have  the  best  quality.  The  second  is  the 
number  of  satellites  in  view.  The  RMS  error  is  lower  when  more  satellites 
are  tracked. 

The  results  of  this  research  show  that  GPS  is  applicable  for  attitude 
determination  in  a  number  of  ways.  The  carrier  phase  approach  allows  for 
its  use  on  spinning  satellites  when  very  good  accuracy  is  required.  For 
missions  with  a  requirement  of  only  coarse  attitude  knowledge,  the  SNR 
technique  is  applicable  and  easily  applied.  It  can  also  be  used  for  initial¬ 
ization  of  more  precise  instruments  or  as  an  emergency  backup.  Its  big¬ 
gest  advantage  is  that  it  does  not  need  a  specialized  attitude  receiver.  If  a 
GPS  navigation  receiver  is  on-board,  then  SNR-based  attitude  estimation 
can  take  place. 

There  are  some  important  areas  where  this  research  can  be 
expanded  in  the  future.  The  phase-based  method  needs  to  be  evaluated  on 
spacecraft  data.  Its  performance  should  be  evaluated  in  terms  of  the  esti¬ 
mation  and  control  requirements  of  an  actual  spinning  satellite  mission. 
Additional  data  is  also  needed  for  the  SNR-based  method.  Measurements 
need  to  be  collected  for  a  two  antenna  configuration  on  a  three-axis  satel- 
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lite.  The  future  SNOE  satellite  mission  might  provide  the  necessary  data 
to  evaluate  this  technique  on  a  spinning  satellite  with  the  antenna 
mounted  on  the  side  of  the  spacecraft.  Finally,  the  performance  of  the 
SNR-based  method  should  be  compared  in  more  detail  with  the  require¬ 
ments  of  various  ambiguity  resolutions  techniques. 
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